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Abstract 

Measurements of the double-differential production cross-section in the range of momentum 
0.5 GeV/c < p < 8.0 GeV/c and angle 0.025 rad < 8 < 0.25 rad in interactions of charged pions on 
beryllium, carbon, aluminium, copper, tin, tantalum and lead are presented. These data represent 
the first experimental campaign to systematically measure forward pion hadroproduction. 

The data were taken with the large acceptance HARP detector in the T9 beam line of the CERN 
PS. Incident particles, impinging on a 5% nuclear interaction length target, were identified by an 
elaborate system of beam detectors. The tracking and identification of the produced particles was 
performed using the forward spectrometer of the HARP detector. Results are obtained for the double- 
differential cross-sections d 2 a/dpdfl mainly at four incident pion beam momenta (3 GeV/c, 5 GeV/c, 
8 GeV/c and 12 GeV/c). The measurements are compared with the GEANT4 and MARS Monte 
Carlo simulation. 
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1 Introduction 



The HARP experiment [TJ at the CERN PS was designed to measure hadron yields from a large range of 
nuclear targets and for incident particle momenta from 1.5 GeV/c to 15 GeV/c. With incident protons, 
this corresponds to a momentum region of great interest for neutrino beams and far from the coverage 
by earlier dedicated hadroproduction experiments [HSU- Measurements with incident charged pions 
are relevant for the calculation of cosmic-ray muon and neutrino fluxes using extended air shower simula- 
tions. These data are also important to simulate re-interactions in particle detectors and neutrino beam 
production targets. At present such calculations rely on models with large uncertainties given the lack 
of experimental data in this energy region. 

Covering an extended range of solid targets in the same experiment, it is also possible to perform system- 
atic comparison of hadron production models with measurements at different incoming beam momenta 
over a large range of target atomic number A. 

This paper presents our final measurements of the double-differential cross-section, d 2 a^ /dpdfl for 
forward production by incident charged pions of 3 GeV/c, 5 GeV/c, 8 GeV/c and 12 GeV/c momentum 
impinging on a thin beryllium, carbon, aluminium, copper, tin, tantalum or lead target of 5% nuclear 
interaction length (A/). In addition, high statistics data at 8.9 GeV/c (12.9 GeV/c) on a thin beryllium 
(aluminium) target are presented. 

Some HARP results on pion production in the forward region with incident protons have already been 
published in papers [6l [3 [8] . 

Pion production data at low momenta and in the forward region from incident pions are extremely 
scarce (SJ. HARP is the first experiment to provide a high statistics data set, taken with many different 
targets, full particle identification and large-acceptance detector. The collected statistics, for the different 
nuclear targets, are reported in Table Q] 

Other data with an incident pion beam, up to higher momenta, has been collected by the Fermilab E970- 
MIPP experiment [TO] and more data at other momenta would be available with its proposed upgrade 

[TJJ. 

The paper is organized as follows. In subsection ll.ll we briefly describe the HARP detector, while section[2] 
is devoted to the main features of the analysis procedure. Section [3] presents results. It is followed by a 
summary presented in section [H 



1.1 Experimental apparatus 

The HARP experiment makes use of a large-acceptance spectrometer consisting of a forward and large- 
angle detection system. The HARP detector is shown in Fig. [TJ A detailed description of the experimental 
apparatus can be found in Ref . |12j . The forward spectrometer - based on five modules of large area drift 
chambers (NDC1-5) [T3J and a dipole magnet complemented by a set of detectors for particle identification 
(PID): a time-of-flight wall (TOFW) [14], a large Cherenkov detector (CHE) and an electromagnetic 
calorimeter (ECAL) - covers polar angles up to 250 mrad. The muon contamination of the beam is 
measured with a muon identifier consisting of thick iron absorbers and scintillation counters. The large- 
angle spectrometer - based on a Time Projection Chamber (TPC) and Resistive Plate Chambers (RPCs) 
located inside a solenoidal magnet - has a large acceptance in momentum and angular range for the pions 
relevant to the production of the muons in a neutrino factory. For the analysis described here only the 
forward spectrometer and the beam instrumentation are used. 

The HARP experiment, located in the T9 beam of the CERN PS, took data in 2001 and 2002. The 
momentum definition of the T9 beam is known with a precision of the order of 1% [15] . 

The target is placed inside the inner field cage (IFC) of the TPC, in an assembly that can be moved in 
and out of the solenoid magnet. The targets used for the measurements reported here have a cylindrical 
shape with a nominal diameter of about 30 mm. Their thickness is equivalent to about 5% Aj. 
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A set of four multi-wire proportional chambers (MWPCs) measures the position and direction of the 
incoming beam particles with an accuracy of ~1 mm in position and ?»0.2 mrad in angle per projection. 
A beam time-of- flight system (BTOF) measures the time difference of particles over a 21.4 m path- 
length. It is made of two identical scintillation hodoscopes, TOFA and TOFB (originally built for the 
NA52 experiment [TB]), which, together with a small target-defining trigger counter, TDS (also used for 
the trigger), provide particle identification at low momenta. This provides separation of pions, kaons 
and protons up to 5 GeV/c and determines the initial time at the interaction vertex (t ). The timing 
resolution of the combined BTOF system is about 70 ps. A system of two N2-filled Cherenkov detectors 
(BCA and BCB) is used to tag electrons at low energies and pions at higher energies. The electron 
and pion tagging efficiency is found to be close to 100%. At the beam energy used for this analysis the 
Cherenkov counters select all particles lighter than protons, while the BTOF is used to reject ions. A set 
of trigger detectors completes the beam instrumentation. 

The selection of beam pions is performed as described in [12] . At 3 GeV/c BCB gas pressure is set to 
tag (BCA is evacuated to reduce multiple scattering of the beam) while the TOFs are capable of 
resolving pions from protons. At 5 GeV/c pions are separated from protons by using joint information 
from the TOFs and one of the Chcrcnkovs (usually BCB), while the gas pressure in the other Cherenkov 
(BCA) is set to tag e^. At higher momenta the contamination drops below 1% and the TOF system 
becomes unable to separate pions from protons efficiently: the beam Cherenkov detectors are used for 
7r-p separation. At the highest beam energy the gas pressure in one of the Cherenkov detectors has been 
set to distinguish protons from kaons. This allows us to estimate the kaon component to be negligible in 
the pion beam selection at all momenta. 



Table 1: Total number of events and tracks used in the various nuclear 5% Ai target data sets and the 
number of incident pions on target as calculated from the pre-scaled incident beam triggers. Numbers 
are for incident ir + (in parenthesis for incident 7r~) in units of 10 3 events. 



Data set (GeV/c) 




3 


5 


8 


8.9 


12 


12.9 


Total DAQ events 


(Be) 


1113 (2233) 


1296 (1798) 


1935 (1585) 


5868 


1207 (1227) 






(C) 


1345 (1831) 


2628 (1279) 


1846 (1399) 




1062 (646) 






(Al) 


1159 (1523) 


1789 (920) 


1707 (1059) 




619 (741) 


4713 




(Cu) 


624 (3325) 


2079 (1805) 


2089 (1615) 




745 (591) 






(Sn) 


1637 (1972) 


2828 (1625) 


2404(1408) 




1803(937) 






(Ta) 


1783 (994) 


2084(1435) 


1965(1505) 




866(961) 






(Pb) 


1911 (1282) 


2111(2074) 


2266 (1496) 




487(1706) 




Accepted beam pions with 


(Be) 


246 (731) 


384 (914) 


341 (826) 


1278 


76 (693) 




forward interaction 


(C) 


257 (299) 


754 (530) 


358 (772) 




41 (352) 






(Al) 


213(486) 


523 (308) 


335 (611) 




27 (435) 


332 




(Cu) 


168 (1185) 


611 (850) 


397 (966) 




33 (347) 






(Sn) 


467 (778) 


819 (732) 


481 (804) 




79 (584) 






(Ta) 


561 (426) 


600 (671) 


388 (893) 




37 (536) 






(Pb) 


611 (473) 


594 (997) 


444 (896) 




20 (839) 




Final state n~ 


(Be) 


.9 (24.9) 


3.5 (32.1) 


6.3 (31.0) 


26.7 


2.0 (27.5) 




selected with PID 


(C) 


.7 (9.1) 


5.9 (15.0) 


5.9 (24.2) 




1.0 (13.7) 






(Al) 


.5 (10.3) 


4.0 (6.8) 


5.6 (13.4) 




.7 (13.4) 


9.8 




(Cu) 


.3 (18.9) 


4.3 (19.0) 


6.0 (25.1) 




.9 (12.9) 






(Sn) 


.7 (9.1) 


4.8 (14.4) 


6.9 (16.3) 




1.9k (19.4) 






(Ta) 


.6 (3.7) 


3.1 (11.5) 


4.9 (16.2) 




.9k (16.6) 






(Pb) 


.6 (3.4) 


3.0 (16.0) 


5.2 (16.3) 




.5k (25.5) 




Final state n + 


(Be) 


8.9 (1.3) 


14 (4.8) 


13.8 (10.6) 


50.0 


3.4 (13.1) 




selected with PID 


(C) 


8.0 (.5) 


23.8 (2.5) 


13.3 (8.0) 




1.7 (6.7) 






(Al) 


5.2 (.5) 


15 (1.3) 


12.1 (5.6) 




1.2 (7.0) 


15.1 




(Cu) 


2.8 (1.6) 


14.8 (3.5) 


12.2 (9.6) 




1.4 (6.6) 






(Sn) 


5.9 (.6) 


16.8 (2.8) 


13.3 (6.1) 




3.1 (10.2) 






(Ta) 


5.3 (.2) 


10.8 (2.2) 


9.4 (6.0) 




1.4 (8.7) 






(Pb) 


5.0 (.2) 


9.8 (2.9) 


10.3 (6.1) 




.7 (13.1) 
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Drift Chambers 




ECAL 



Figure 1: Schematic layout of the HARP detector. The convention for the coordinate system is shown 
in the lower-right corner. 

A downstream trigger in the forward scintillator trigger plane (FTP) was required to record the event, 
accepting only tracks with a trajectory outside the central hole (60 mm) which allows beam particles to 
pass. The trigger counter covers the spectrometer acceptance fully and has a high efficiency (>99.8%). 

The length of the accelerator spill was 400 ms with a typical intensity of 15 000 beam particles per spill. 
The average number of events recorded by the data acquisition ranged from 300 to 350 per spill. 

The absolute normalization of the number of incident pions was performed using 'incident-particle' trig- 
gers. These are triggers where the same selection on the beam particle was applied but no selection on 
the interaction was performed. The rate of this trigger was down-scaled by a factor 64. 



2 Data Analysis 

2.1 Event and particle selection 

A detailed description of the experimental techniques used for data analysis in the HARP forward spec- 
trometer can be found in Ref. [5J 117) . Further details of the improved analysis techniques can be found 
in [5j [7J. For the current analysis we have used identical reconstruction and PID algorithms, while at 
the final stage of the analysis the unfolding technique introduced as UFO in [5J has been applied. This 
technique has already been described in details in Ref. [TH] . 

At the first stage of the analysis a beam pion is selected using the beam time of flight system (TOF-A, 
TOF-B) and the Cherenkov counters (BCA, BCB) as described in section [Ol We always require time 
measurements in TOF-A, TOF-B and/or TDS to be present which are needed for calculating the arrival 
time of the beam pion at the target. 

Secondary track selection criteria, described in [Jj, are optimized to ensure the quality of momentum 
reconstruction and a clean time-of-flight measurement while maintaining a high reconstruction efficiency. 

The background induced by interactions of beam particles in the materials outside the target is measured 
by taking data without a target in the target holder ("empty target data"). These data are subject to 
the same event and track selection criteria as the standard data sets. 
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To take into account this background the number of particles of the observed type (n + , n~ ) in the "empty 
target data" are subtracted bin-by-bin (momentum and angular bins) from the number of particles of 
the same type. The uncertainty induced by this method is discussed in section 12.31 and labeled "empty 
target subtraction" . The event statistics is summarized in Table [TJ 

The negative beam consists of e~ and 7r _ (with a dominant fraction of 7r~), while the positive beam 
is dominated by protons at high momenta and by tt + at low momenta 0. The kaon background is 
estimated to be <0.5% and is neglected in the analysis, while the muon contamination (tagged by the 
muon identifier) is around 3% and is subtracted. 

The small fraction of pions in the positively charged beam explains the significantly different statistics 
of the data-sets with incident tt~ and 7r + . The relatively small fraction of accepted final state pions 
compared to the accepted beam pions is due to the stringent quality cuts applied in the analysis. The 
major loss is due to the small acceptance of the dipole magnet in the vertical plane, and the fact that 
only focused particles are accepted for analyses. The latter requirement reduces the statistics by a factor 
two but improves the quality of reconstruction and reduces the background [6]. 



2.2 Cross-section calculation 



The cross-section is calculated as follows 



p'9'..a' 



where 



• dpdn (pi , Oj ) is the cross-section in mb/(GeV/c sr) for the particle type a (p, n + or 7r ) for each 
true momentum and angle bin (pi,6j) covered in this analysis; 

• N a (jo'i , 9j ) is the number of particles of type a' in bins of reconstructed momentum p\ and angle 
9'j in the raw data; 

• ■Mp'eapie'a' i s the correction matrix which accounts for efficiency and resolution of the detector; 

• ]vjpt ' Iv^T an d Ap An are normalization factors, namely: 

N ^P t is the number of target nuclei per unit aread; 

A/pot is the number of incident beam particles on target (particles on target); 
Api and Aflj are the bin sizes in momentum and solid angle, respectively 0. 



We do not make a correction for the attenuation of the incoming beam in the target, so that strictly 
speaking the cross-sections are valid for Ai = 5% targets. 

The calculation of the correction matrix M™n _,_/«/„/ is a rather difficult task. Various techniques are 

Pi " j " j ^ 

described in the literature to obtain this matrix. As discussed in Ref. [5] for the p-Al analysis of HARP 
data at 12.9 GeV/c, two complementary analyses have been performed to cross-check internal consistency 
and possible biases in the respective procedures. A comparison of both analyses shows that the results 
are consistent within the overall systematic error [5] . 

In the first method - called "Atlantic" in [5] - the correction matrix M^. ap , e i a , is decomposed into 
distinct independent contributions, which are computed mostly using the data themselves. The second 

1 the proton fraction in the incoming beam goes from 35% at 3 GeV/c to about 92% at 12 GeV/c 
2 A - atomic mass, Na - Avogadro number, p - target density and t - target thickness 
3 Api = p™ ax - pf ll \ AQj = 27r(cos(6» min ) - cos(6» max )) 
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method - called "UFO" in [5] - is the unfolding method introduced by D'Agostini [18 0. This method 
has been used in the recent HARP publications [5J 030 [201 1211 H3 and it is also applied in the analysis 
described here (see [7] for additional information). 

The Monte Carlo simulation of the HARP setup is based on the GEANT4 package [231 HI- The detec- 
tor materials are accurately described in this simulation as well as the relevant features of the detector 
response and the digitization process. All relevant physics processes are considered, including multiple 
scattering, energy loss, absorption and re-interactions. The simulation is independent of the beam mo- 
mentum and particle type because it only generates for each event exactly one secondary particle of a 
specific particle type inside the target material and propagates it through the complete detector. A small 
difference (at the few percent level) is observed between the efficiency calculated for events simulated 
with the single-particle Monte Carlo and with a simulation using a multi-particle hadron-production 
model. A similar difference is seen between the single-particle Monte Carlo and the efficiencies measured 
directly from the data. The single-particle Monte Carlo predicts an efficiency more in agreement with 
the efficiency measured with the data. A momentum-dependent correction factor determined using the 
efficiency measured with the data is applied to take this into account. The track reconstruction algo- 
rithms used in this analysis and the simulation are identical to the ones used for the 7r + production in 
p-Be collisions [6]. A detailed description of the corrections and their magnitude can be found there. 

The reconstruction efficiency (inside the geometrical acceptance) is larger than 95% above 1.5 GeV/c and 
drops to 80% at 0.5 GeV/c. The requirement of a match with a TOFW hit has an efficiency between 
90% and 95% above 1.0 GeV/c and is responsible for the drop in efficiency towards lower momenta. 
The electron veto rejects about 1% of the pions and protons below 3 GeV/c with a remaining electron 
background of less than 0.5%. Below Cherenkov threshold the TOFW separates pions and protons with 
negligible background and an efficiency of «98% for pions. Above Cherenkov threshold the efficiency for 
pions is greater than 99% with only 1.5% of the protons mis-identified as a pion. The kaon background 
in the pion spectra is smaller than 1%. 

The absorption and decay of particles is simulated by the Monte Carlo. The generated single particle 
can re-interact and produce background particles by hadronic or electromagnetic processes, thus giving 
rise to tracks in the forward spectrometer ( "tertiaries" ) . In such cases also the additional measurements 
are entered into the migration matrix thereby taking into account the combined effect of the generated 
particle and any secondaries it creates. The absorption correction is on average 20%, approximately 
independent of momentum. Uncertainties in the absorption of secondaries in the dipole spectrometer 
material are taken into account by a variation of 10% of this effect in the simulation. The effect of pion 
decay is treated in the same way as the absorption and is 20% at 500 MeV/c and negligible at 3 GeV/c. 

The uncertainty in the production of background due to tertiary particles is larger. The average correction 
is ~10% and up to 20% at 1 GeV/c. The correction includes re-interactions in the detector material as 
well as a small component coming from re-interactions in the target. The validity of the generators used 
in the simulation was checked by an analysis of HARP data with incoming protons and charged pions 
on aluminium and carbon targets at lower momenta (3 GeV/c and 5 GeV/c). A 30% variation of the 
secondary production was applied. The average empty-target subtraction amounts to «20%. 

Owing to the redundancy of the tracking system downstream of the target the reconstruction efficiency is 
very robust under the usual variations of the detector performance during the long data taking periods. 
Since the momentum is reconstructed without making use of the upstream drift chamber module (which 
is more sensitive in its performance to the beam intensity) the reconstruction efficiency is uniquely 
determined by the downstream system. No variation of the overall efficiency has been observed. The 
performance of the TOFW and CHE system have been monitored to be constant for the data taking 
periods used in this analysis. The calibration of the detectors was performed on a day-by-day basis. 

4 The unfolding method tries to put in correspondence the vector of measured observables (such as particle momentum, 
polar angle and particle type) a: mea s with the vector of true values xtrue using a migration matrix: x meas = M m i gr X xtrue- 
The goal of the method is to compute a transformation (correction matrix) to obtain the expected values for xtrue from 
the measured ones. The most simple and obvious solution, based on simple matrix inversion M~} , is usually unstable and 
is dominated by large variances and strong negative correlations between neighbouring bins. In the method of D' Agostini, 
the correction matrix M UFO tries to connect the measurement space (effects) with the space of the true values (causes) 
using an iterative Bayesian approach, based on Monte Carlo simulations to estimate the probability for a given effect to be 
produced by a certain cause. 



2.3 Error estimation 



Different types of sources induce systematic errors for the analysis described here: track yield corrections 
(~ 5%), particle identification (~ 0.1%), momentum and angular reconstruction (~ 0.5%) 0. The domi- 
nant error is due to track yield correction, in particular due to the subtraction of tertiary particles, with 
smaller contributions from particle absorption, empty target subtraction and reconstruction efficiency. 
The strategy to calculate these systematic errors and the different methods used for their evaluation 
are described in [7J where the numerical values are tabulated. An additional source of error is due to 
misidentified secondary kaons, which are not considered in the particle identification method used for this 
analysis and are subtracted on the basis of a Monte Carlo simulation, as in [7J. No explicit correction 
is made for pions coming from decays of other particles created in the target, as they give a very small 
contribution according to the selection criteria applied in the analysis. 

Examples of experimental uncertainties are shown in Figure [5] for Be target with incident ir~ for ir^ 
production at two incident beam momenta (5 and 8 GeV/c). They are very similar for the other beam 
momenta and targets. Going from lighter (Be, C) to heavier targets (Ta, Pb) the corrections for 7r° 
conversion and absorption/tertiaries increase slightly. 

The overall normalization has an uncertainty of ~ 2% and is mainly due to the uncertainty in the efficiency 
that beam pions counted in the normalization actually hit the target, with smaller components from the 
target density and the beam particle counting procedure. On average the total integrated systematic 
error is around 5 — 6%, with a differential bin-to-bin systematic error of the order of 10 — 11%, to be 
compared with a statistical integrated (bin-to-bin differential) error of ~ 2 — 3% (~ 10 — 13%). Systematic 
and statistical errors are roughly of the same order. 



3 Results 

The measured double-differential cross-sections for the production of tt + and 7r _ in the laboratory system 
as a function of the momentum and the polar angle for each incident beam momentum are shown in Fig- 
ures [3] to [9] for targets from Be to Pb. The error bars shown are the square-roots of the diagonal elements 
in the covariance matrix, where statistical and systematic uncertainties are combined in quadrature. The 
correlation of the statistical errors (introduced by the unfolding procedure) are typically smaller than 
20% for adjacent momentum bins and even smaller for adjacent angular bins. The correlations of the 
systematic errors are larger, typically 80% for adjacent bins. The overall scale error (< 2%) is not shown. 
The results of this analysis are also tabulated in Appendix A. We have not included the scale errors in the 
tables to make it possible to calculate e.g. integrated particle ratios taking the scale errors into account 
only when applicable, i.e. when different beams are compared. 

The dependence of the averaged pion yields on the atomic number A is shown from Fig. [10] to Fig. 1131 
The pion yields, averaged over two angular regions (0.05 rad < 8 < 0.15 rad and 0.15 rad < 9 < 0.25 rad) 
and four momentum regions (0.5 GeV/c < p < 1.5 GeV/c, 1.5 GeV/c < p < 2.5 GeV/c, 2.5 GeV/c < p < 
3.5 GeV/c and 3.5 GeV/c < p < 4.5 GeV/c), are shown for four different beam momenta. One observes 
a smooth behaviour of the averaged yields. The most striking feature is the difference between same-sign 
and opposite-sign pion production compared to the beam particle. Here one observes a stronger beam 
energy dependence for opposite-sign production. The A-dependence is similar for n~ and 7r + production 
if one compares opposite-sign and same-sign production, although one observes that especially at lower 
beam momenta the opposite-sign pion production saturates earlier towards higher A. 

In the following we will show only some comparisons of the HARP data with publicly available Monte 
Carlo simulations: GEANT4 [23] (version 9.1p02) and MARS [25] . not previously tuned to our data-sets, 
by using different models. The comparison will be shown for a limited set of plots and only for the Be 
(Figures [T4l to [T5]) and Ta (Figures \W\ to |2"3"]) targets, as examples of a light and a heavy target. 

At intermediate energies (up to 5-10 GeV), GEANT4 uses two types of intra- nuclear cascade models: the 

5 The quoted error in parenthesis refers to fractional error of the integrated cross-section in the kinematic range covered 
by the HARP experiment 
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Figure 2: Systematic errors as a function of momentum of the outgoing pions for the particular case 
of 5 and 8 GeV/c tt~ interacting on a beryllium target. Plots in the upper left panel are for incident 
5 GeV/c 7r~ producing 7r~, plots in the upper right panel are for 5 GeV/c n~ producing 7r + , plots in 
the lower left panels are for 8 GeV/c ir~ producing tt~ and plots in the lower right panel are for 8 
GeV/c 7T _ producing tt + . Total systematic error (grey solid line) are displayed together with the main 
components: black short-dashed line for absorption+tertiaries interactions, black dotted line for track 
efficiency and target pointing efficiency, black dot-dashed line for 7r° subtraction, black solid line for 
momentum scale+resolution and angle scale, grey short-dashed line for PID. 



Bertini model [27l [28] (valid up to ~ 10 GeV) and the binary model [29] (with a validity range up to ~ 3 
GeV). Both models treat the target nucleus in detail, taking into account density variations and tracking 
in the nuclear field. The binary model is based on hadron collisions with nucleons, giving resonances that 
decay according to their quantum numbers. The Bertini model is based on the cascade code reported 
in [30] and hadron collisions are assumed to proceed according to free-space partial cross-sections and 
final-state distributions measured for the incident particle types. 

At higher energies, instead, two parton string models, the quark-gluon string (QGS) model [37] [31] and 
the Fritiof (FTP) model [31] are used, in addition to a High Energy Parametrized model (HEP) derived 
from the high energy part of the GHEISHA code used inside GEANT3 [32] . 

The parametrized models of GEANT4 (HEP and LEP) are intended to be fast, but conserve energy and 
momentum on average and not event by event. 

A realistic GEANT4 simulation is built by combining models and physics processes into what is called a 
"physics list". In high energy calorimetry the two most commonly used are the QGSP physics list, based 
on the QGS model, the pre-compound nucleus model and some of the Low Energy Parametrized (LEP) 
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model [^| and the LHEP physics list based on the parametrized LEP model and HEP models. 

Currently is also popular the QGSP-BERT physics list in which the Bertini model is used at 3 and 5 
GeV/c and the QGSP model at 8 and 12 GeV/c. 



Also this model, at low energy, has its root in the Gheisha code inside GEANT3. 

11 



HARP 7T~ — Be tt + 



HARP it' -Be 7T 



1000- 

100: 

10 



^ 1 

^ 1000|r 

a 
ft 

^ 100 

b 

w 

10 

1 



50-100 mrad 



o 12 
• 8 

a 5 
▼ 3 



150-200 mrad 



f ft 



100-150 mrad 

: 1000 
: ~ 100 

; o 
> 

; ^ j > 1 

200-250 mrad- d 

: — 1000 

: a 

■ -a 

- a 

*- : ^ 100 



b 

CM 
T3 



— X- 



0246802468 
P [GeV/c] 
HARP 7T + -Be tt + 



10 
1 



50-100 mrad 


100-150 mrad 










o 12 
• 8 

a 5 
▼ 3 




150-200 mrad 


200-250 mrad 











0246802468 
P [GeV/c] 
HARP 7T + -Be 7T" 





100-150 mrad- 




0246802468 
P [GeV/c] 



0246802468 
P [GeV/c] 



Figure 3: Be (top) and 7r + — Be (bottom) differential cross-sections for different incoming beam 

momenta (from 3 to 12 GeV/c). Left panels show the tt + production, while right panels show the n" 
production. In the top right corner of each plot the covered angular range is shown in mrad. 
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Figure 4: 7r~— C (top) and 7r + — C (bottom) differential cross-sections for different incoming beam mo- 
menta (from 3 to 12 GeV/c). Left panels show the w + production, while right panels show the ir~ 
production. In the top right corner of each plot the covered angular range is shown in mrad. 
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Figure 5: 7r~— Al (top) and 7r + — Al (bottom) differential cross-sections for different incoming beam mo- 
menta (from 3 to 12 GeV/c). Left panels show the w + production, while right panels show the ir~ 
production. In the top right corner of each plot the covered angular range is shown in mrad. 
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Figure 6: 7r~ — Cu (top) and ir + — Cu (bottom) differential cross-sections for different incoming beam 
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Figure 7: 7r~— Sn (top) and 7r+— Sn (bottom) differential cross-sections for different incoming beam 
momenta (from 3 to 12 GeV/c). Left panels show the n + production, while right panels show the tt~ 
production. In the top right corner of each plot the covered angular range is shown in mrad. 
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Figure 8: 7r~— Ta (top) and 7r + — Ta (bottom) differential cross-sections for different incoming beam 
momenta (from 3 to 12 GeV/c). Left panels show the ir + production, while right panels show the tt~ 
production. In the top right corner of each plot the covered angular range is shown in mrad. 
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Figure 9: 7r~— Pb (top) and n + — Pb (bottom) differential cross-sections for different incoming beam 
momenta (from 3 to 12 GeV/c). Left panels show the ir + production, while right panels show the tt~ 
production. In the top right corner of each plot the covered angular range is shown in mrad. 
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Figure 10: The dependence on the atomic number A of the 7r~ production yields in 7r~-Be, 7r _ -C, 7r~- 
Al, 7r~-Cu, 7r~-Sn, 7r~-Ta, 7r~-Pb interactions averaged over two forward angular region (0.05 rad < 
9 < 0.15 rad and 0.15 rad < 9 < 0.25 rad) and four momentum regions (0.5 GeV/c < p < 1.5 GeV/c, 
1.5 GeV/c < p < 2.5 GeV/c, 2.5 GeV/c < p < 3.5 GeV/c and 3.5 GeV/c < p < 4.5 GeV/c), for the four 
different incoming beam momenta (from 3 GeV/c to 12 GeV/c). 
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Figure 11: The dependence on the atomic number A of the tt + production yields in 7r~-Be, ir~-C, ir~- 
Al, 7r~-Cu, 7r~-Sn, 7r _ -Ta, 7r~-Pb interactions averaged over two forward angular region (0.05 rad < 
9 < 0.15 rad and 0.15 rad < 9 < 0.25 rad) and four momentum regions (0.5 GcV/c < p < 1.5 GeV/c, 
1.5 GcV/c < p < 2.5 GeV/c, 2.5 GcV/c < p < 3.5 GeV/c and 3.5 GcV/c < p < 4.5 GeV/c), for the four 
different incoming beam momenta (from 3 GeV/c to 12 GeV/c). 
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Figure 12: The dependence on the atomic number A of the tt~ production yields in tt + Be, ir + - C, ir + - 
Al, 7r+-Cu, 7r + -Sn, 7r + -Ta, 7r + -Pb interactions averaged over two forward angular region (0.05 rad < 
9 < 0.15 rad and 0.15 rad < 6 < 0.25 rad) and four momentum regions (0.5 GcV/c < p < 1.5 GeV/c, 
1.5 GeV/c < p < 2.5 GeV/c, 2.5 GcV/c < p < 3.5 GeV/c and 3.5 GcV/c < p < 4.5 GeV/c), for the four 
different incoming beam momenta (from 3 GcV/c to 12 GeV/c). 
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Figure 13: The dependence on the atomic number A of the ir + production yields in 7r + -Be, 7r + -C, 7r + 
Al, 7r + -Cu, 7r + --Sn, 7r + -Ta, 7r+-Pb interactions averaged over two forward angular region (0.05 rad < 
9 < 0.15 rad and 0.15 rad < 6 < 0.25 rad) and four momentum regions (0.5 GcV/c < p < 1.5 GeV/c, 
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different incoming beam momenta (from 3 GeV/c to 12 GeV/c). 
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Figure 14: Comparison of HARP double-differential tt^ 1 cross-sections for 7r-Be at 3 GeV/c with GEANT4 
and MARS MC predictions, using several generator models. 



23 



y lit~Be->n~X 50-100 mrad ~| 



>10' 
a 
O 

m 



ST - I ' • 



■D 

a. 

"5 





|;i~Be-»7[~X 100-150 mrad | 



i .... i .... i .... j_ 



|7i'Be->7['X 150-200 mrad | 




i .... i .... i .... j_ 



1 7i~ Be -» 7i" X 200-250 mrad | 



I 1 11 1 I ■ 1 ■ 1 I ■ 1 1 ■ I J I 



10 3 



7rBe->K*X 50-100 mrad | 



1 7i- Be -> Ti* X 100-1 50 mrad | 



|ti- Be -> 71* X 1 50-200 mrad | 



| it- Be -» 7i* X 200-250 mrad | 




[-•-HH 









10 3 



10 2 



10 



r Be^i'X 50-100 mrad | 



1 7i* Be -» 7i' X 1 00-1 50 mrad | 



|7i+Be^7TX 150-200 mrad | 



| 7t* Be -> ti" X 200-250 mrad | 









i .... i .... i .... i 




i .... i .... i . . i 



10 3 



10 2 



10 



ti* Be -> ti* X 50-100 mrad 



7i*Be->7t*X 100-150 mrad 



7t*Be^7t*X 150-200 mrad 



ti* Be -> ti* X 200-250 mrad 




■ i i ■ I i i i ■ I i i ■ ■ I i ■ ■ ■ I i ■ 



I .... I .... I .... I 



I .... I .... I 



I .... I .... I .... I 



12 3 4 



12 3 4 

p [GeV/c] 



Figure 15: Comparison of HARP double-differential tt^ 1 cross-sections for 7r-Be at 5 GeV/c with GEANT4 
and MARS MC predictions, using several generator models. 
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Figure 16: Comparison of HARP double-differential 7r ± cross-sections for n -Be at 8 GeV/c with 
GEANT4 and MARS MC predictions, using several generator models. 
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Figure 17: Comparison of HARP double-differential ^ cross-sections for 7r + -Be at 8 GeV/c with 
GEANT4 and MARS MC predictions, using several generator models. 
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Figure 18: Comparison of HARP double-differential n ± cross-sections for n-Be at 12 GeV/c with 
GEANT4 and MARS MC predictions, using several generator models. 



27 



The MARS code system |25] uses as basic model an inclusive approach to multiparticle production 
originated by R. Feynman in 1974. At each interaction vertex, a particle cascade tree is constructed by 
using a fixed number of particles with weights that in the simplest case are equal to the partial mean 
multiplicity of that event. Above 3 GeV phenomenological particle production models are used, while 
below 5 GeV a cascade-exciton model [53] combined with the Fermi break-up model, the coalescence 
model, an evaporation model and a multifragmentation extension are used instead. 

HARP measurements have been compared to the available models, first renormalizing MC predictions to 
data and then computing the x 2 between models and data themselves. Normalization factors are shown 
in Tables [2] and [3J An absolute comparison yields x 2 values that are too large to be meaningful. Over 
the full energy range covered by the HARP experiment, the best comparison is obtained with the MARS 
Monte Carlo which is on average 10% below the data, with a maximum shift of 30%. This may be owing 
to the fact that MARS is using different models in different energy regions, equivalent to use a collection 
of models as implemented in the "physics lists" of GEANT4. Unfortunately, all models have very bad 
X 2 compared with the data, showing significant inconsistencies in the shape distribution. This is worse 
than the situation with incident protons and may be explained by the fact that in the latter case some 
experimental data were available for tuning of the Monte Carlo simulations. This was not the case with 
incident pions in the momentum range studied here (below 15 GeV/c). 

At higher energies the FTP model (from GEANT4) and the MARS model describe the data better, while 
at the lowest energies the Bertini and binary cascade models (from GEANT4) seem more appropriate. 
Parametrized models, such as LHEP from GEANT4, show important discrepancies. 



Table 2: Normalization factors data over MC for some models, incident n + . 
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Table 3: Normalization factors data over MC for some models, incident ir 
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The full set of HARP data, taken with targets spanning the full periodic table of elements, with small 
total errors and full coverage of the solid angle in a single detector will definitely help the validation 
of models used in hadronic simulations in the difficult energy range between 3 GeV/c and 15 GeV/c of 
incident momentum with incident pions, where data were previously missing. 
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Figure 19: Comparison of HARP double-differential tt^ 1 cross-sections for 7r-Ta at 3 GeV/c with GEANT4 
and MARS MC predictions, using several generator models. 
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Figure 20: Comparison of HARP double-differential tt^ 1 cross-sections for 7r-Ta at 5 GeV/c with GEANT4 
and MARS MC predictions, using several generator models. 
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Figure 21: Comparison of HARP double-differential ^ cross-sections for 7r -Ta at 8 GeV/c with 
GEANT4 and MARS MC predictions, using several generator models. 
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Figure 22: Comparison of HARP double-differential ir* 1 cross-sections for 7r + -Ta at 8 GeV/c with 
GEANT4 and MARS MC predictions, using several generator models. 
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Figure 23: Comparison of HARP double-differential n ± cross-sections for 7r-Ta at 12 GeV/c with 
GEANT4 and MARS MC predictions, using several generator models. 
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4 Summary and conclusions 



Double-differential cross-sections for the production of positive and negative pions in the kinematic range 
0.5 GeV/c < < 8 GcV/c and 0.025 rad < 6 n < 0.25 rad from the interactions of 3, 5, 8 and 12 GeV/c 
7T^~ on beryllium, carbon, aluminium, copper, tin, tantalum and lead targets of 5% interaction length 
thickness have been reported. 

We should stress that the HARP data are the first precision measurements with incoming charged pion 
beams in this kinematic region and may have a major impact on the tuning of Monte Carlo generators at 
low energies. Quantitative comparisons with a variety of generators currently in use have been presented. 

The pion yield averaged over different momentum and angular ranges increase smoothly with the atomic 
number A of the target and with the energy of the incoming pion beam. 



5 Acknowledgments 

We gratefully acknowledge the help and support of the PS beam staff and of the numerous technical 
collaborators who contributed to the detector design, construction, commissioning and operation. In 
particular, we would like to thank G. Barichello, R. Brocard, K. Burin, V. Carassiti, F. Chignoli, D. Con- 
venti, G. Decreuse, M. Delattre, C. Detraz, A. Domeniconi, M. Dwuznik, F. Evangelisti, B. Friend, A. Ia- 
ciofano, I. Krasin, D. Lacroix, J.-C. Lcgrand, M. Lobello, M. Lollo, J. Loquct, F. Marinilli, R. Mazza, 
J. Mulon, L. Musa, R. Nicholson, A. Pepato, P. Petev, X. Pons, I. Rusinov, M. Scandurra, E. Uscnko, 
R. van der Vlugt, for their support in the construction of the detector and P. Dini for his contribu- 
tion to Monte Carlo production. The collaboration acknowledges the major contributions and advice 
of M. Baldo-Ceolin, L. Linssen, M.T. Muciaccia and A. Pullia during the construction of the experi- 
ment. The collaboration is indebted to V. Ableev, F. Bergsma, P. Binko, E. Boter, M. Calvi, C. Cavion, 
M.Chizov, A. Chukanov, A. DeSanto, A. DeMin, M. Doucet, D. Dullmann, V. Ermilova, W. Flegel, 
Y. Hayato, A. Ichikawa, O. Klimov, T. Kobayashi, D. Kustov, M. Laveder, M. Mass, H. Mcinhard, 
A. Menegolli, T. Nakaya, K. Nishikawa, M. Paganoni, F. Paleari, M. Pasquali, M. Placentino, V. Ser- 
diouk, S. Simonc, P.J. Solcr, S. Troquereau, S. Ueda, A. Valassi and R. Veenhof for their contributions 
to the experiment. 

We acknowledge the contributions of V. Ammosov, G. Chelkov, D. Dedovich, F. Dydak, M. Gostkin, 
A. Guskov, D. Khartchenko, V. Koreshev, Z. Kroumchtein, I. Nefcdov, A. Semak, J. Wotschack, V. Zaets 
and A. Zhemchugov to the work described in this paper. 

The experiment was made possible by grants from the Institut Interuniversitaire des Sciences Nucleair- 
es and the Interuniversitair Instituut voor Kernwetcnschappen (Belgium), Ministerio de Educacion y 
Ciencia, Grant FPA2003-06921-c02-02 and Generalitat Valenciana, grant GV00-054-1, CERN (Geneva, 
Switzerland), the German Bundesministerium fur Bildung und Forschung (Germany), the Istituto Nazio- 
nale di Fisica Nucleare (Italy), INR RAS (Moscow), the Russian Foundation for Basic Research (grant 
08-02-00018), the Particle Physics and Astronomy Research Council (UK) and the Swiss National Science 
Foundation, in the framework of the SCOPES programme. We gratefully acknowledge their support. 



References 

[1] M. G. Catanesi et at, HARP Collaboration, "Proposal to study hadron production for the neutrino 
factory and for the atmospheric neutrino flux", CERN-SPSC/99-35 (1999). 

[2] G. Ambrosini et al. [NA56 Collaboration], Eur. Phys. J. C 10 (1999) 605; 
G. Ambrosini et al. [NA56 Collaboration], Phys. Lett. B 420 (1998) 225; 
G. Ambrosini et al. [NA56 Collaboration], Phys. Lett. B 425 (1998) 208. 

[3] HW. Athcrton, CERN 80-07, August 1980. 



34 



M. Bonesini and A. Guglielmi, Phys. Rep. 433 (2006) 66. 

M. G. Catanesi et al. [HARP Collaboration], "Measurement of the production cross-section of 
positive pions in p-Al collisions at 12.9 GeV/c", Nucl. Phys. B 732 (2006) 1 [arXiv:hep-ex/0510039| . 

M. G. Catanesi et al. [HARP Collaboration], "Measurement of the production cross-section of 
positive pions in the collision of 8.9 GeV/c protons on beryllium", Eur. Phys. J. C 52 (2007) 29 
[arXiv:hep ex/0702024| . 

M. G. Catanesi et al. [HARP Collaboration], "Measurement of the production cross-sections of tt^ 
in p-C and ti^-C interactions at 12 GeV/c" , Astroparticle Physics 29/4 (2008) 257, larXiv:0802.0657l 
[astro-ph] . 

M. G. Catanesi et al. [HARP Collaboration], " Forward production in p-C-2 and P-N2 interactions 
at 12 GeV/c", Astroparticle Physics 30/3 (2008) 124. larXiv:0807.1025l [hep-ex] . 

H.R.Blieden et al, Phys. ReV. Dll, 14 (1975). 

MIPP Collaboration, |http://ppd.fnal.gov/experiments/290771 

R. Raja, "The Main Injector particle production Experiment (MIPP) at Fermilab", Jour. Phys. 
Conf. Series 9 (2005) 303. 

D. Isenhower et al, Proposal to update the MIPP experiment, Fermilab Proposal P-960. 

M. G. Catanesi et al. [HARP Collaboration], "The HARP Detector at the CERN PS", Nucl. 
Instrum. Meth. A 571 (2007) 527. 

M. Anfreville et al, "The drift chambers of the NOMAD experiment", Nucl. Instrum. Meth. A 481 
(2002) 339 |arXiv:hep-ex/0104012| . 



M. Baldo-Ceolin et al, "The Time-Of-Flight TOFW Detector Of The HARP Experiment: Con- 
struction And Performance" , Nucl. Instrum. Meth. A 532 (2004) 548; 

M. Bonesini et al., "Construction of a Fast-laser Calibration System for the HARP TOF counter 
wall", IEEE Trans. Nucl. Science 50 (2003) 541. 

L. Durieu, A. Mueller and M. Martini, PAC-2001-TPAH142 Presented at IEEE Particle Accelerator 
Conference (PAC2001), Chicago, Illinois, 18-22 Jun 2001; 
L. Durieu et al, Proceedings of PAC'97, Vancouver, (1997); 
L. Durieu, O. Fernando, CERN PS/PA Note 96-38. 

K. Pretzl et al., Invited talk at the "International Symposium on Strangeness and Quark Matter", 
Crete, (1999) 230. 

M. G. Catanesi et al. [HARP Collaboration], "Particle identification algorithms for the HARP 
forward spectrometer", Nucl. Instrum. Meth. A 572 (2007) 899. 

G. D'Agostini, DESY 94-099, ISSN 0418-9833, 1994. 
G. D'Agostini, Nucl. Instrum. Meth. A 362 (1995) 487. 

M. G. Catanesi et al., [HARP Collaboration], "Measurement of the Production of charged Pions 
by Protons on a Tantalum Target", Eur. Phys. J. C 51 (2007) 787. larXiv:0706.1600l [hep-ex] 

M. G. Catanesi et al., [HARP Collaboration], "Large-angle production of charged pions by 3 GeV/c- 
12 GeV/c protons on carbon, copper and tin targets" , Eur. Phys. C 53 (2008) 177. larXiv:0709.3464l 
[hep-ex] , 

M. G. Catanesi et al., [HARP Collaboration], "Large-angle production of charged pions by 3 GeV/c- 
12.9 GeV/c protons on beryllium, aluminium and lead targets", Eur. Phys. C 54 (2008) 37, 
larXiv:0709.3458l [hep-ex] . 

M. G. Catanesi et al., [HARP Collaboration], "Large-angle production of charged pions with inci- 
dent protons on nuclear targets as measured in the HARP experiment", Phys. ReV. C77 (2008) 
055207. 



35 



[23] S. Agostinelli et al. [GEANT4 Collaboration], Nucl. Instrum. Meth. A 506 (2003) 250. 

[24] V. Ivanchenko, "Simulation and recent results of the HARP experiment", NSS Conference Record, 
2006, IEEE, Vol3, pp. 1660-1665 

[25] N.V. Mokhov, S.I. Striganov, "MARS overview", FERMILAB-CONF-07-008-AD, 2007. 

[26] D.H. Wright et al, AIP Conf. Proc. 867 (2006) 479. 

[27] D.H. Wright et al, AIP Conf. Proc. 896 (2007) 11. 



[28] A. Heikkinen et al. e-print physics/0306008 



[29] G. Folger, V. Ivanchenko and H.P. Wellisch, Eur. Phys. J. A21 (2004) 407. 
[30] H.W. Bertini, P. Guthrie, Nucl. Phys. A169 (1971). 
[31] G. Folger and H.P. Wellisch, e-print |physics/0306007| 
[32] H. Fesefeld, Technical report PITHA 85-02, Aachen, 1985. 
[33] S.G. Mashnik et al, LANL report LA-UR-05-7321, 2005. 



36 



A Cross-section data 



The following tables report the measured differential cross-section for positive and negative pion pro- 
duction in interactions of 3, 5, 8 and 12 GeV/c momentum charged pions on different types of nuclear 
targets. The data are presented in the kinematic range of 0.5 GeV/c < < 8 GeV/c and 0.05 rad 
<6^< 0.25 rad. The overall normalization uncertainty ( < 2%) is not included in the reported errors. 

Results at higher incoming beam momenta, from 8 GeV/c to 12.9 GcV/c, are also presented extending 
to a lower value of the polar angle 9 (0.025 rad), using a finer binning. This second set of tables contains 
additional data for collisions of 12.9 GeV/c 7r ± on Al and 8.9 GeV/c tt^ 1 on Be. 

The cross-section values in some bins of the following tables are omitted and replaced by the symbol 
* ± *, due to instabilities in the unfolding procedure coming from the low statistics available. 
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Table 4: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 I '(dpdSl) , for 7r~ -Be interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariance matrix are given. 
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Table 5: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa n /(dpdCt), for n~-Be interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariance matrix are given. 
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Table 6: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 a 77 j (dpdVl) , for 7r + -Be interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, #min — ^ < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 7: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d 2 I '(dpdSl) , for 7r + -Be interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < "max) bin, where p and 8 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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0.001±0.004 


0.002±0.006 
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Table 8: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 a n / (dpdQ.) , for ir~-C interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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± 


0.002±0.000 
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0.09±0.02 
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0.20±0.03 






1.00 
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0.14±0.02 
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0.15±0.03 






2.00 
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0.010±0.006 


0.039±0.008 


0.069±0.012 
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2.50 
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0.018±0.005 
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Table 9: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa n /(dpdCt), for n~-C interactions at 3,5,8,12 GeV/c. Each row refers to a different {p m [ n < p < 
Pmax, #min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 10: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 a- 71 " / (dpdCt) , for w + -C interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, #min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix arc given. 
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0.119±0.014 


0.18±0.02 


0.18±0.06 






1.50 


2.00 


0.12±0.02 


0.105±0.011 


0.17±0.02 


0.13±0.05 






2.00 


2.50 


0.09±0.02 


0.080±0.010 


0.13±0.02 


0.12±0.07 






2.50 


3.00 


± 


0.055±0.007 


0.080±0.013 


0.10±0.08 






3.00 


3.50 


± 


0.042±0.010 


0.076±0.012 


0.11±0.04 






3.50 


4.00 


± 


0.031±0.007 


0.043±0.008 


0.06±0.03 






4.00 


5.00 


± 


0.022±0.012 


0.025±0.005 


0.022±0.015 
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0.002±0.004 
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3.00 


± 


0.034±0.007 


0.055±0.014 


0.05±0.05 






3.00 


3.50 


± 


0.016±0.004 


0.040±0.011 


0.016±0.016 
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4.00 
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0.012±0.005 


0.032±0.009 


0.010±0.030 
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0.008±0.004 


0.023±0.008 
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5.00 
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± 


± 
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Table 11: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7T+-C interactions at 3,5,8,12 GeV/c. Each row refers to a different {p m [ n < p < 
Pmax, ^min — ^ < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 12: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 a- 71 " / (dpdfl) , for 7T~-A1 interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m [ n < p < 
Pmax, #min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 13: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7r~-Al interactions at 3,5,8,12 GeV/c. Each row refers to a different {p m [ n < p < 
Pmax, ^min — ^ < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 14: HARP results for the double-differential 7r + production cross-section in the laboratory system, 
d 2 a- 71 " / (dpdfl) , for 7r + -Al interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m [ n < p < 
Pmax, #min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 15: HARP results for the double-differential it~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7T+-A1 interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m j n < p < 
Pmax, ^min — ^ < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 16: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 a- 71 " / (dpdVl) , for 7r~-Cu interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m ; rl < p < 
Pmax, #min — ^ < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 17: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7r~-Cu interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m j n < p < 
Pmax, ^min — ^ < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariance matrix are given. 
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9 ^n 
z.ou 


n 99_i_n n/i 
U.zzztU.U4 


0.31zb0.03 


0.63±0.05 


U.ozztU.U ( 






9 ^n 
z.ou 


q nn 
o.uu 


zt 


0.31zb0.03 


0.40±0.03 


n fiQ-i-n n^ 

U.OoztU.UO 






q nn 
o.uu 


Q c:n 

o.ou 


1 

zt 


0.34zb0.03 


0.31±0.03 


n /io_i_n c\k 
U.4yztU.U0 






q p;n 
o.ou 


a nn 

4.UU 


_l_ 
zt 


0.21zb0.03 


0.18±0.03 


n QQ_i_n nQ 
U.ooztU.Uo 






/i nn 

4.UU 


^ nn 

O.UU 


_l_ 
zt 


0.21zb0.04 


0.118±0.014 


n 9^-i-n nQ 

U.ZOztU.Uo 






£ nn 
O.UU 


o.ou 


_l_ 

zt 


± 


0.059±0.007 


n 1 nn-i-n ni q 
U.lUUztU.Ulo 






a p;n 
o.ou 


s nn 

o.UU 


1 

zt 


± 


± 


U.UoOztU.UUu 


U. 10 


n on 
u.zu 


n p;n 

U.OU 


i nn 
1 .uu 


n 99-i-n c\a 

U.ZZztU.U4 


0.42zb0.05 


0.88±0.09 


i n^-i-n i q 

l.UOztU.lo 






i nn 

l.UU 


i t^n 
1 .ou 


n 91 -i~n nQ 
U.zlztU.Uo 


0.33±0.04 


0.76±0.07 


n QA-i-n no 
u.ooztu.uy 






i ^n 

l.OU 


9 nn 

Z.UU 


n 1 tx_|_n 09 
U.lOztU.UZ 


0.27±0.03 


0.51±0.05 


U.OOztU.U ( 






9 nn 
z.UU 


9 ^n 
z.ou 


n i 7-i-n no 
U.l / ztU.Uz 


0.19zb0.02 


0.35±0.03 


U.OUztU.UO 






9 p;n 
z.ou 


q nn 
o.uu 


i 

zt 


0.12±0.02 


0.20±0.02 


n Q9_i_n c\a 
U.ozztU.U4 






q nn 
o.uu 


q ^n 
o.ou 


zt 


0.073±0.012 


0.114=b0.015 


n i Q-i-n nQ 
u.iyztu.uo 






q p;n 
o.OU 


a nn 
4.UU 


_l_ 
zt 


0.039±0.008 


0.097=b0.012 


n i c;-i-n n9 
U.lOztU.UZ 






4.00 


5.00 


± 


0.036±0.007 


0.049=b0.007 








5.00 


6.50 


zb 


± 


0.019=b0.003 


0.037±0.007 






6.50 


8.00 


zb 


± 


± 


0.009±0.003 


0.20 


0.25 


0.50 


1.00 


0.17=b0.03 


0.39±0.06 


0.59=b0.07 


0.72±0.10 






1.00 


1.50 


0.11=b0.02 


0.18±0.03 


0.39±0.05 


0.38±0.06 






1.50 


2.00 


0.11zb0.02 


0.14±0.03 


0.31=b0.05 


0.42±0.07 






2.00 


2.50 


0.056=b0.011 


0.07±0.02 


0.35=b0.05 


0.31±0.06 






2.50 


3.00 


zb 


0.047±0.011 


0.23=b0.04 


0.16±0.04 






3.00 


3.50 


zb 


0.040±0.009 


0.13=b0.02 


0.13±0.03 






3.50 


4.00 


zb 


0.032±0.007 


0.070=b0.013 


0.11±0.02 






4.00 


5.00 


zb 


0.026±0.006 


0.041=b0.008 


0.08±0.02 






5.00 


6.50 


zb 


± 


0.013=b0.004 


0.029±0.010 






6.50 


8.00 


zb 


± 


zb 


0.006±0.006 
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Table 18: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 a- 71 " / (dpdTt) , for 7r + -Cu interactions at 3,5,8, 12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, #min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 



mm 
(rad) 


^ITiaX 

(rad) 


ft 

^min 

(GeV/c) 


(GeV/c) 


3 GeV/c 


d 2 a n+ /{dpdQ) 
(barn/(GcV/c sr)) 
5 GeV/c 8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.16±0.06 


a ooj_n nc 


A C 0_l_A AA 


0.74±0.25 






1.00 


1.50 


0.08±0.02 


a ooj_n c\a 
U.oz±U.U4 


a £in_Ln ac 


0.81±0.21 






1.50 


2.00 


0.15±0.02 


a qt_i_a no 


A C7_I_A ac 


0.75±0.20 






2.00 


2.50 


0.43±0.05 


a o a _i_ a no 


A KC_I_A ac 


0.68±0.18 






2.50 


3.00 


± 


a /in_LA c\a 
U.4y±U.U4 


A K A _I_A A C 


0.77±0.19 






3.00 


3.50 


± 


a ca_i_a ac 


A A 0_l_A A E 


0.68±0.15 






3.50 


4.00 


± 


a 7Q i n c\a 
U. ( o±U.U4 


A /1C_I_A AK 


0.71±0.16 






4.00 


5.00 


± 


a coin oo 
U.Oo±U. zz 


A /I Q_I_A AO 


0.55±0.10 






5.00 


6.50 


± 


i 

± 


A A 1 _1_A AC 

U.41±U.Uo 


0.41±0.07 






6.50 


8.00 


± 


i 

± 


1 

± 


0.18±0.05 


0.100 


0.150 


0.50 


1.00 


0.12±0.05 


a a K_i_n ac 


A C /I _I_A nc 


0.86±0.27 






1.00 


1.50 


0.18±0.05 


a oq_i_a nc 


A C 1 _J_ A AC 


0.95±0.25 






1.50 


2.00 


0.28±0.05 


a /1a_i_a nc 
U.4U±U.Uo 


A CQ_I_A f17 


0.79±0.20 






2.00 


2.50 


0.51±0.07 


a o/i_i_n c\a 
U.o4zhU.U4 


A EC_I_A AC 


0.56±0.16 






2.50 


3.00 


± 


n o a _i_n a a 
U.o4±U.U4 


A O 1 _l_ A A A 


0.65±0.16 






3.00 


3.50 


± 


A OO 1 n AO 


A O 1 _l_ A A A 


0.42±0.13 






3.50 


4.00 


± 


A OT_l_A AO 

U.z /±U.Uo 


A 01 _I_A AO 


0.23±0.07 






4.00 


5.00 


± 


a 97+n 1 


A 1 CIA (Y\ 


0.23±0.07 






5.00 


6.50 


± 


± 


0.066±0.014 


0.07±0.03 






6.50 


8.00 


± 


± 


± 


0.03±0.02 


0.150 


0.200 


0.50 


1.00 


0.21±0.07 


0.50±0.07 


0.74±0.10 


0.85±0.30 






1.00 


1.50 


0.22±0.05 


0.32±0.04 


0.55±0.06 


0.58±0.19 






1.50 


2.00 


0.23±0.05 


0.26±0.03 


0.39±0.05 


0.68±0.19 






2.00 


2.50 


0.14±0.04 


0.19±0.03 


0.30±0.04 


0.72±0.21 






2.50 


3.00 


± 


0.13±0.02 


0.24±0.04 


0.28±0.11 






3.00 


3.50 


± 


0.09±0.02 


0.15±0.03 


0.23±0.11 






3.50 


4.00 


± 


0.07±0.02 


0.069±0.015 


0.16±0.08 






4.00 


5.00 


± 


0.05±0.03 


0.043±0.010 


0.14±0.07 






5.00 


6.50 


± 


± 


0.017±0.007 


0.021±0.024 






6.50 


8.00 


± 


± 


± 


0.007±0.010 


0.200 


0.250 


0.50 


1.00 


0.18±0.09 


0.28±0.05 


0.53±0.08 


0.41±0.22 






1.00 


1.50 


0.07±0.03 


0.23±0.04 


0.33±0.06 


0.34±0.26 






1.50 


2.00 


0.10±0.04 


0.14±0.03 


0.26±0.05 


0.21±0.17 






2.00 


2.50 


0.05±0.03 


0.15±0.03 


0.25±0.05 


0.14±0.20 






2.50 


3.00 


± 


0.09±0.02 


0.14±0.04 


0.124±0.253 






3.00 


3.50 


± 


0.041±0.012 


0.10±0.03 


0.037±0.087 






3.50 


4.00 


± 


0.021±0.009 


0.06±0.02 


0.034±0.108 






4.00 


5.00 


± 


0.021±0.011 


0.038±0.014 


0.037±0.124 






5.00 


6.50 


± 


± 


0.018±0.011 


0.012±0.036 






6.50 


8.00 


± 


± 


± 


0.003±0.102 
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Table 19: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7r+-Cu interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m ; n < p < 
Pmax, ^min — ^ < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 



ft ■ 




^min 


pmax 




d 2 a 7r+ /(dpdCl) 




[TtXil j 


yidid J 




frvv /A 

V / C) 




(barn/(GeV/c sr)) 












O ^jG v / c 


5 GcV/c 


8 GeV/c 


1Z KjQ V / C 


n n^n 
U.UOU 


n i nn 


n c^n 
U.OU 


i nn 
l.UU 


n 11 _i_n n/i 
U.llztU.U4 


0.27±0.04 


0.51±0.07 


n QQ-i_n i c\ 

U.OOztU.lu 






i nn 


1 .OU 


n i Q-i~n n/i 

U.10ztU.U4 


0.24zb0.03 


0.54±0.06 


n 79-i-n 91 

U. ( ZztU.Zl 






l.OU 


9 nn 

Z.UU 


n i Q-i~n n^ 

U.lOztU.UO 


0.20zb0.03 


0.41±0.05 


U.DoztU.iy 






o nn 
z.UU 


9 r^n 
z.ou 


nil -i~n n/i 
U.llztU.U4 


0.19zb0.02 


0.38±0.04 


n a i-\-f\ i a 

U.4 ( ztU.14 






9 p;n 
z.ou 


q nn 
o.uu 


zt 


0.16zb0.02 


0.27±0.03 


n s^-i-n 9n 

U.oOztU.ZU 






q nn 
o.uu 


q f^n 
o.ou 


_i_ 
zt 


0.16zb0.02 


0.25±0.03 


n 90-i-n no 
u.zyztu.uy 






q ^n 
o.ou 


/i nn 


_i_ 
zt 


0.18zb0.02 


0.18±0.02 


n a i -i~n 1 9 

U.41ztU.lZ 






/i nn 
4.UU 


x. nn 
O.UU 


_l_ 
zt 


0.07zb0.02 


0.17±0.02 


n 90-i-n ns 
u.zyztu.uo 






£ nn 

o.uu 


c;n 
D.OU 


_l_ 
zt 


± 


0.065±0.008 


n i 9-i~n c\a 
U.lzztU.U4 






D.OU 


s nn 
o.uu 


1 

zt 


± 


± 


n no-i~n c\a 
u.uyztu.u4 


n i nn 


U.10U 


n c^n 
U.OU 


i nn 
l.UU 


n 97-i-n c\7 

U.Z / ztU.U / 


0.45zb0.06 


0.65±0.09 


i i n~i~n q^ 

1 .lUztU.OO 






i nn 

l.UU 


i ^n 

1 .OU 


n i/i _j_n n/i 

U.14ztU.U4 


0.23zb0.03 


0.46±0.05 


n Qc;-i-n 99 

U.oOztU. zz 






l.OU 


9 nn 

Z.UU 


n i n~i~n nQ 
U.lUztU.Uo 


0.21zb0.03 


0.38±0.04 


n qi -i~n i 9 
U.olztU.lz 






9 nn 

Z.UU 


9 ^n 
z.ou 


n n/i -i~n n9 
U.U4ztU.Uz 


0.18zb0.02 


0.29±0.04 


n a Q-i~n i ^ 
U.4oztU.10 






9 p;n 
z.ou 


q nn 
o.uu 


zt 


0.100±0.015 


0.24±0.03 


n /i/i -un i a 

U.44ztU.14 






q nn 
o.uu 


q f^n 
o.ou 


1 

zt 


0.077±0.013 


0.18±0.02 


n Q/1 _i_n 1 1 
U.o4ztU.ll 






q p;n 
o.ou 


/i nn 

4.UU 


_l_ 
zt 


0.031±0.006 


0.082±0.014 


n QA-i-fi i 9 
U.ouztU.lz 






/i nn 
4.UU 


x. nn 
o.uu 


_l_ 
zt 


0.022±0.007 


0.073±0.012 


n i Q-i~n c\« 
U.loztU.Uu 






£ nn 
O.UU 


^ c;n 
D.OU 


_l_ 
zt 


± 


0.018±0.004 


n n/i -i~n nQ 
U.U4ztU.Uo 






p;n 

D.OU 


s nn 
o.UU 


1 

zt 


± 


± 


n nn9-i~n nnQ 
U.UUzztU.UUo 




n onn 
u.zuu 


n p;n 

U.OU 


i nn 

1 .UU 


nil -i~n c\a 

U.llZtU.U4 


0.42zb0.06 


0.58±0.08 


n 70-i-n 97 
u. < yztu.z ( 






i nn 

l.UU 


i p;n 

1 .OU 


n i n~i~n n/i 
U.lUztU.U4 


0.23zb0.03 


0.42±0.06 


n C\(\-\-C\ 9Q 

U.DuztU. zo 






i ^n 

l.OU 


9 nn 

Z.UU 


nil -un n/i 

U.llztU.U4 


0.23±0.03 


0.33±0.04 


n qi 4~n 1 q 

U.OlztU.lO 






9 nn 
z.UU 


9 r^n 
z.ou 


n ni Q-i~n ni q 
U.UlyztU.Ulo 


0.11±0.02 


0.28±0.04 


n 9Q_i_n 1 1 
U.zOztU.ll 






9 p;n 
z.ou 


q nn 
o.uu 


i 

zt 


0.064±0.013 


0.16±0.03 


n i c\-\-(\ no 
u.iuztu.uy 






q nn 
o.uu 


q ^n 
o.ou 


_l_ 
zt 


0.01d±0.004 


0.08±0.02 


n 99-i-n i n 

U.ZZztU.lU 






q c^n 
o.OU 


/i nn 
4.UU 


_l_ 
zt 


0.012±0.005 


0.037±0.010 


n i Q-i~n no 
u.ioztu.uy 






4.00 


5.00 


-±- 


0.005±0.003 


0.020±0.006 


o 1 1 +n 07 






5.00 


6.50 


zb 


± 


0.004±0.002 


0.03±0.03 






6.50 


8.00 


zb 


± 


± 


± 


0.200 


0.250 


0.50 


1.00 


0.17=b0.07 


0.33±0.06 


0.60±0.09 


0.58±0.26 






1.00 


1.50 


0.06=b0.03 


0.20±0.04 


0.51±0.08 


0.44±0.23 






1.50 


2.00 


0.08=b0.04 


0.09±0.02 


0.29±0.05 


0.72±0.37 






2.00 


2.50 


0.009=b0.011 


0.024±0.007 


0.13±0.03 


0.16±0.18 






2.50 


3.00 


zb 


0.014±0.005 


0.043±0.012 


0.07±0.11 






3.00 


3.50 


zb 


0.017±0.007 


0.020±0.009 


0.050±0.065 






3.50 


4.00 


zb 


0.008±0.004 


0.008±0.006 


0.010±0.031 






4.00 


5.00 


zb 


0.005±0.004 


0.005±0.007 


0.019±0.109 






5.00 


6.50 


zb 


± 


0.002±0.004 


0.016±0.090 






6.50 


8.00 


zb 


± 


± 


0.003±0.068 
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Table 20: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 a- 71 " / (dpdCt) , for 7T~-Sn interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, #min — ^ < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 



mm 


^max 


^mm 






d 2 a^ 


/(dpdQ.) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c sr)) 












3 GeV/c 


O Vjii V / c 


O \jfki V / C 


12 GeV/c 


0.05 


0.10 


0.50 


1.00 


0.11±0.03 


U.ozzhU.Uo 


n (^Q-i_n os 
U.DoztU.Uo 


1.15±0.12 






1.00 


1.50 


0.07±0.02 


U.o4ztU.U4 


n QQ_i_n 07 


1.25±0.09 






1.50 


2.00 


0.14±0.03 


o oo-i-o ri/i 


n t'O-lo o^ 


1.02±0.08 






2.00 


2.50 


0.11±0.03 


o i Q-i-H oo 


n c^o-i-o o^ 


0.93±0.07 






2.50 


3.00 


± 


o i 0*3 




0.75±0.05 






3.00 


3.50 


± 


o 1 7-4-0 no 

U.l / rtU.UZ 


n ^fi-i-o o^ 

U.oOztU.UO 


0.69±0.06 






3.50 


4.00 


± 


n 1 7-i-H oo 


n oi 4-o O/i 


0.50±0.04 






4.00 


5.00 


± 


U.UOztU.Uo 


n i p;-i-0 09 


0.36±0.03 






5.00 


6.50 


± 


zt 


o nso-i-0 oi q 


0.19±0.02 






6.50 


8.00 


± 


1 

zt 


i 

zt 


0.081±0.011 


0.10 


0.15 


0.50 


1.00 


0.18±0.03 


U.4ZztU.UO 


n Q/1 -LO 1 


1.21±0.12 






1.00 


1.50 


0.07±0.02 


U.oO=tU.U4 


n «o_i_n 07 


1.23±0.11 






1.50 


2.00 


0.08±0.02 


o oQ-i-H oq 


n fii -i-o c\« 


1.07±0.09 






2.00 


2.50 


0.12±0.03 


U.ZlztU.Uo 


n c^o-i-o o^ 

U.OUztU.UO 


0.87±0.08 






2.50 


3.00 


± 


n i a -t-n no 


n Q/i-i-o O/i 


0.65±0.06 






3.00 


3.50 


± 


n oo-i-o no 
U.UyzhU.Uz 


O OO-l-O OQ 


0.55±0.05 






3.50 


4.00 


± 


n n^fi^n nno 


n i 74-0 OQ 


0.41±0.04 






4.00 


5.00 


± 


n ooo^o nno 
u.uzzztu.uuy 


O 07-1-0 OO 

U.U / ±U.Uz 


0.21±0.02 






5.00 


6.50 


± 


± 




0.067±0.011 






6.50 


8.00 


± 


± 


± 


0.008±0.003 


0.15 


0.20 


0.50 


1.00 


0.22±0.05 


0.45±0.07 


0.92±0.11 


1.22±0.13 






1.00 


1.50 


0.10±0.03 


0.24±0.03 


0.57±0.06 


0.86±0.08 






1.50 


2.00 


0.040±0.012 


0.24±0.03 


0.50±0.05 


0.82±0.07 






2.00 


2.50 


0.05±0.02 


0.14±0.03 


0.34±0.04 


0.52±0.06 






2.50 


3.00 


± 


0.06±0.02 


0.19±0.03 


0.33±0.04 






3.00 


3.50 


± 


0.031±0.010 


0.12±0.02 


0.20±0.03 






3.50 


4.00 


± 


0.021±0.008 


0.07±0.02 


0.09±0.02 






4.00 


5.00 


± 


0.004±0.004 


0.027±0.007 


0.053±0.010 






5.00 


6.50 


± 


± 


0.005±0.002 


0.013±0.004 
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Table 21: HARP results for the double-differential it~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7r~-Sn interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 22: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 a- 71 " / (dpdCt) , for 7r + -Sn interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, #min — ^ < "max) bin, where p and 8 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 23: HARP results for the double-differential 7r~ production cross-section in the laboratory system, 
d 2 a n /(dpdCl), for 7r + -Sn interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < "max) bin, where p and 8 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 24: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 cr 71 " / (dpdTt) , for 7r~-Ta interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 



mm 
(rad) 


(rad) 


^mm 
(GeV/c) 


pITiaX 

(GeV/c) 


3 GeV/c 


d 2 a 7r+ /(dpdQ) 
(barn/(GcV/c sr)) 
5 GeV/c 8 GeV/c 


12 GeV/c 


0.05 


0.10 


0.50 


1.00 


0.08±0.03 


n /io_i_n a*s 
U.42±U.Uo 


1 A£!_I_A 1 1 

l.UOzhU.ll 


1.60±0.17 






1.00 


1.50 


0.06±0.03 


a oo_i_a nc 


1 OT_l_A AA 

L.zf ±u.uy 


1.67±0.13 






1.50 


2.00 


0.13±0.05 


a oa_i_a c\a 


1 A A _l_ A AT 

l.U4±U.U ( 


1.27±0.09 






2.00 


2.50 


0.17±0.06 


a o a _i_ a a a 
U.z4±U.U4 


A 7fl_Ln f\T 

U. fU±U.U ( 
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3.00 


± 
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3.50 


± 
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0.89±0.07 
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± 
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U.l t rtU.Uo 
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± 
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0.43±0.04 
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± 


i 

± 
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± 


i 

± 


1 

± 
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0.15 


0.50 
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± 
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0.78±0.07 
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± 
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± 
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6.50 


8.00 


± 


± 


± 


0.014±0.003 


0.15 


0.20 


0.50 


1.00 


0.11±0.05 


0.45±0.08 


1.22±0.14 


1.78±0.19 






1.00 


1.50 


0.12±0.05 


0.29±0.05 


0.64±0.07 


0.98±0.10 






1.50 


2.00 


0.06±0.02 


0.24±0.04 


0.51±0.06 


0.90±0.09 






2.00 


2.50 


0.06±0.04 


0.18±0.04 


0.41±0.05 


0.71±0.08 






2.50 


3.00 


± 


0.04±0.02 


0.27±0.04 


0.45±0.06 






3.00 


3.50 


± 


0.040±0.015 


0.16±0.03 


0.27±0.04 






3.50 


4.00 


± 


0.009±0.007 


0.10±0.02 


0.14±0.03 






4.00 


5.00 


± 


0.001±0.002 


0.024±0.007 


0.046±0.012 






5.00 


6.50 


± 


± 


0.004±0.002 


0.011±0.005 






6.50 


8.00 


± 


± 


± 


0.002±0.001 


0.20 


0.25 


0.50 


1.00 


0.25±0.11 


0.53±0.10 


1.01±0.14 


1.24±0.17 






1.00 
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0.15±0.07 


0.51±0.09 
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1.50 


2.00 


0.07±0.05 


0.30±0.06 


0.79±0.12 
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2.00 
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0.02±0.03 


0.20±0.04 


0.42±0.08 
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0.05±0.02 


0.22±0.05 
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3.00 
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± 
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5.00 


6.50 


± 


± 
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0.001±0.001 






6.50 


8.00 


± 


± 


± 


0.000±0.013 
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Table 25: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7r~-Ta interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < #max) bin, where p and are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariance matrix are given. 
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O.UU 


zt 


0.31zb0.05 


0.21±0.03 
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U.OozbU.U4 
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O.UU 


o.ou 
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zb 


± 
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zb 


± 


± 
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U.OU 
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1.31±0.12 


1 QO_i_n i 7 
l.oyzbU.l ( 






l.OU 


9 nn 

Z.UU 


n 99-i-n n7 

U.ZZztU.U ( 


0.39±0.06 


0.88±0.09 


1. lUzbU. 1 1 






9 nn 
z.UU 


9 ^n 
z.ou 


n 1 A-i~n n/i 
U.luztU.U4 


0.26zb0.05 


0.53±0.06 


n si -i~n ns 
U.olzbU.Uo 






9 p;n 
z.ou 


^ nn 
o.uu 


zt 


0.12±0.02 


0.31±0.04 


U.4uzbU.UO 






q nn 
o.uu 


q ^n 
o.ou 


zt 


0.11±0.02 


0.17±0.03 


n Af\-\-f\ n^ 

U.4UzbU.UO 






q p;n 
o.OU 


a nn 
4.UU 


zb 


0.07±0.02 


0.16±0.02 


n oi -J_n no 
U.zlzbU.Uo 






4.00 


5.00 


± 


0.06±0.02 


0.079±0.014 


16+0 02 






5.00 


6.50 


zb 


± 


0.034±0.007 


0.047±0.010 






6.50 


8.00 


± 


± 


± 


0.012±0.004 


0.20 


0.25 


0.50 


1.00 


0.17=b0.07 


0.50±0.09 


1.10±0.14 


1.22±0.16 






1.00 


1.50 


0.11=b0.05 


0.30±0.05 


0.40±0.06 


0.58±0.08 






1.50 


2.00 


0.06=b0.03 


0.23±0.05 


0.39±0.08 


0.64±0.10 






2.00 


2.50 


0.03=b0.02 


0.15±0.04 


0.28±0.05 


0.49±0.08 






2.50 


3.00 


± 


0.08±0.02 


0.21±0.04 


0.32±0.07 






3.00 


3.50 


± 


0.10±0.03 


0.15±0.03 


0.21±0.04 






3.50 


4.00 


zb 


0.050±0.014 


0.08±0.02 


0.14±0.03 






4.00 


5.00 


± 


0.026±0.009 


0.06±0.02 


0.11±0.03 






5.00 


6.50 


± 


± 


0.025±0.010 


0.035±0.010 






6.50 


8.00 


zb 


± 


± 


0.012±0.008 
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Table 26: HARP results for the double-differential 7r + production cross-section in the laboratory system, 
d 2 cr 71 " / (dpdVt) , for 7r + -Ta interactions at 3,5,8,12 GeV/c. Each row refers to a different {p m i n < p < 
Pmax, ^min — ^ < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 



mm 
(rad) 


^ITiaX 

(rad) 


ft 

^min 
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3 GeV/c 
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0.26±0.07 
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0.56±0.18 






4.00 


5.00 


± 


n Q«+n 14 


n OK_i_n 


0.52±0.15 






5.00 


6.50 


± 


± 


0.08±0.02 


0.18±0.07 






6.50 


8.00 


± 


± 


± 


0.10±0.05 


0.150 


0.200 


0.50 


1.00 


0.23±0.07 


0.79±0.13 


1.38±0.20 


1.11±0.49 






1.00 


1.50 


0.28±0.06 


0.42±0.08 


0.70±0.10 


0.78±0.32 






1.50 


2.00 


0.23±0.04 


0.36±0.06 


0.56±0.08 


1.33±0.41 






2.00 


2.50 


0.21±0.04 


0.38±0.07 


0.51±0.08 


0.75±0.30 






2.50 


3.00 


± 


0.14±0.03 


0.28±0.06 


0.22±0.13 






3.00 


3.50 
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0.16±0.04 


0.25±0.04 


0.27±0.16 






3.50 


4.00 


± 


0.12±0.03 


0.21±0.04 


0.15±0.09 






4.00 


5.00 


± 


0.07±0.04 


0.11±0.03 


0.11±0.08 
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0.029±0.012 
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6.50 


8.00 


± 
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± 


0.042±0.045 


0.200 


0.250 


0.50 
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0.30±0.13 


0.44±0.09 


0.69±0.13 


0.74±0.37 






1.00 


1.50 


0.15±0.04 


0.24±0.06 


0.50±0.11 


0.96±0.56 
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2.00 


0.13±0.04 


0.23±0.06 


0.37±0.09 


0.45±0.28 






2.00 


2.50 


0.024±0.014 


0.17±0.05 


0.39±0.09 


1.09±0.53 
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0.09±0.03 


0.23±0.08 


0.48±0.38 






3.00 
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0.02±0.10 
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± 


± 
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0.02±0.20 
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Table 27: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7r+-Ta interactions at 3,5,8,12 GeV/c. Each row refers to a different {p m [ n < p < 
Pmax, ^min — ^ < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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± 
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± 


± 
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U.UzztU.Uz 




u.zuu 
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i nn 
1 .uu 
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0.65±0.11 


1.23±0.18 


1 .DOztU.O ( 






i nn 

l.UU 


i p;n 
1 .ou 


n i fi-i-n n/i 
U.luztU.U4 


0.41zb0.06 


0.75±0.11 


1 91 -i~n QO 
l.zlztU.oy 






i ^n 

l.OU 


9 nn 

Z.UU 


n nG-i-n n°. 
u.uyztu.uo 


0.29±0.05 


0.60±0.09 


1 .lOztU.O ( 






9 nn 
z.UU 


z.ou 


n n/i -i~n n9 
U.U4ztU.Uz 


0.19±0.03 


0.49±0.08 
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z.ou 
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o.uu 
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zt 
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n /i 9-i~n 91 
U.4zztu.zl 






°. nn 
o.uu 
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U.loztU.lz 
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5.00 


± 


0.006±0.004 


0.024±0.010 


1 5+0 1 1 






5.00 


6.50 


zb 


± 


0.003±0.003 


0.017±0.037 






6.50 


8.00 


zb 


± 


± 


0.000±0.001 


0.200 


0.250 


0.50 


1.00 


0.26=b0.08 


0.52±0.11 


1.07±0.18 


1.24±0.54 






1.00 


1.50 


0.18=b0.06 


0.27±0.06 


0.72±0.13 


0.87±0.37 






1.50 


2.00 


0.05=b0.03 


0.11±0.04 


0.35±0.07 


1.61±0.74 






2.00 


2.50 


0.009=b0.008 


0.07±0.02 


0.12±0.04 


1.14±0.60 






2.50 


3.00 


zb 


0.04±0.02 


0.026±0.012 


0.20±0.20 






3.00 


3.50 


zb 


0.008±0.006 


0.016±0.013 


0.12±0.21 






3.50 


4.00 


zb 


O.OOiO.OO 


0.00±0.01 


0.12±0.27 






4.00 


5.00 


zb 


O.OOiO.OO 


0.00±0.01 


0.04±0.19 






5.00 


6.50 


zb 


± 


0.00±0.01 


0.00±0.04 






6.50 


8.00 


zb 


± 


± 


0.00±0.10 
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Table 28: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 a- 71 " / (dpdVl) , for 7r~-Pb interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, #min — ^ < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 



mm 
(rad) 


^max 
(rad) 


^mm 
(GeV/c) 


(GeV/c) 


3 GeV/c 


d 2 a 7r+ /(dpdQ) 
(barn/(GeV/c sr)) 
5 GeV/c 8 GeV/c 


12 GeV/c 


0.05 


0.10 


0.50 


1.00 


0.04±0.02 


n on _l n n£! 


n nn i n i i 

u.yy±u.n 


1.64±0.16 






1.00 


1.50 


0.04±0.02 


n a o_i_n n c 
U.4z±U.Uo 
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1.4b±U.lU 


1.75±0.12 
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U.DO±U.U/ 
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3.00 


± 
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U.41±U.Uo 


0.94±0.07 
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± 


n i n_i_n no 
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U.4Z±U.Uo 


0.96±0.07 






3.50 


4.00 


± 


n i K_i_n no 


n oi _i_n 
U.ol±U.UO 


0.64±0.07 






4.00 


5.00 


± 


n n£!_i_n no 
U.UOztU.Uz 
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0.37±0.03 
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6.50 


± 


i 

± 
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U.lU±U.Uz 


0.20±0.02 






6.50 
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± 


i 

± 


i 

± 


0.074±0.011 


0.10 


0.15 


0.50 


1.00 


0.18±0.06 


U.4o±U.U ( 


1.40±U.l0 


1.89±0.18 






1.00 


1.50 


0.05±0.02 


n on_i_n n c 


n hk _i_n nn 

u.yo±u.uy 


1.59±0.13 






1.50 


2.00 


0.09±0.05 


n o7_lh c\a 
U.zY ±U.U4 


n 7K _l n no 


1.21±0.10 






2.00 


2.50 


0.06±0.03 


n o a _i_n n a 
U.z4±U.U4 


n fiKi n n7 


1.00±0.09 






2.50 


3.00 


± 


n i £!_i_n no 


n o7 i n nc 


0.72±0.07 






3.00 


3.50 


± 


nil _i_n no 
U.ll±U.Uz 


n ok _i_n n a 
U.ZO±U.U4 


0.57±0.06 






3.50 


4.00 


± 


n nn_i_n no 

u.uy±u.uz 


n i 7j_n no 
U.l ( ±U.Uo 


0.45±0.04 






4.00 


5.00 


± 


n m Q-LC] ni o 


n 09 

U.UOZLL/.UZi 


0.24±0.03 






5.00 


6.50 


± 


± 


0.016±0.007 


0.085±0.014 






6.50 


8.00 


± 


± 


± 


0.018±0.004 


0.15 


0.20 


0.50 


1.00 


0.02±0.02 


0.48±0.08 


1.29±0.15 


1.86±0.19 






1.00 


1.50 


0.14±0.07 


0.26±0.04 


0.80±0.08 


0.99±0.09 






1.50 


2.00 


0.020±0.011 


0.22±0.03 


0.68±0.07 


0.94±0.09 






2.00 


2.50 


0.04±0.04 


0.14±0.03 


0.45±0.06 


0.73±0.07 






2.50 


3.00 


± 


0.07±0.02 


0.31±0.05 


0.43±0.05 






3.00 


3.50 


± 


0.07±0.02 


0.15±0.03 


0.37±0.04 






3.50 


4.00 


± 


0.011±0.009 


0.07±0.02 


0.17±0.03 






4.00 


5.00 


± 


0.001±0.001 


0.042±0.010 


0.08±0.02 






5.00 


6.50 


± 


± 


0.008±0.003 


0.012±0.005 






6.50 


8.00 


± 


± 


± 


0.002±0.001 


0.20 


0.25 


0.50 


1.00 


0.20±0.11 


0.63±0.10 


1.01±0.14 


1.44±0.18 






1.00 


1.50 


0.10±0.06 


0.54±0.08 


0.97±0.13 


1.53±0.18 
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2.00 


0.12±0.07 


0.33±0.06 


0.75±0.12 


1.20±0.16 
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2.50 
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0.18±0.04 
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3.00 


± 


0.05±0.02 


0.21±0.05 


0.32±0.05 
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± 


0.016±0.011 
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± 
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± 
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0.024±0.012 
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6.50 


± 


± 


0.002±0.002 


0.003±0.002 






6.50 


8.00 


± 


± 


± 


0.001±0.009 
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Table 29: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7r~-Pb interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < "max) bin, where p and 8 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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3.00 


zb 
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3.00 
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0.14±0.03 
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3.50 


4.00 


zb 
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zb 


0.031±0.010 
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Table 30: HARP results for the double-differential 7r + production cross-section in the laboratory system, 
d?a w /(dpdCl), for 7r+-Pb interactions at 3,5,8,12 GeV/c. Each row refers to a different {p m i n < p < 
Pmax, ^min — " < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariance matrix are given. 





min 


#max 




Pmax 




a a 


/(dpdQ.) 




(rad) 


(rad) 


(GeV/c) 


(GeV/c) 




(barn/(GeV/c sr)) 












3 GeV/c 


5 GcV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.20±0.06 


0.64±0.12 


0.90±0.16 


2.89±1.10 






1.00 


1.50 


0.042±0.011 


0.38±0.07 


0.78±0.10 


0.82±0.42 






1.50 


2.00 


0.093±0.011 


0.43±0.07 


0.77±0.10 


1.54±0.65 






2.00 


2.50 


0.29±0.03 


0.72±0.08 


0.82±0.11 


0.73±0.36 






2.50 


3.00 


± 


0.75±0.07 


0.74±0.10 


1.64±0.53 






3.00 


3.50 


± 


0.89±0.09 


0.61±0.08 


1.55±0.50 






3.50 


4.00 


± 


1.04±0.07 


0.68±0.09 


1.38±0.45 






4.00 


5.00 


± 


0.78±0.26 


0.89±0.07 


0.73±0.23 






5.00 


6.50 


± 


± 


0.72±0.11 








6.50 


8.00 


± 


± 


± 


0.14±0.15 


0.100 


0.150 


0.50 


1.00 


0.08±0.04 


0.87±0.14 


1.40±0.19 


2.46±0.91 






1.00 


1.50 


0.18±0.06 


0.36±0.07 


1.00±0.15 


2.03±0.73 






1.50 


2.00 


0.52±0.08 


0.55±0.08 


1.07±0.21 


0.54±0.31 






2.00 


2.50 


0.69±0.09 


0.52±0.07 


0.89±0.36 


1.20±0.65 






2.50 


3.00 


± 


0.49±0.06 


0.49±0.43 


0.88±0.53 






3.00 


3.50 


± 


0.66±0.08 


0.57±0.12 


0.70±0.37 






3.50 


4.00 


± 


0.57±0.06 


0.40±0.06 


0.53±0.26 






4.00 


5.00 


± 


0.35±0.13 


0.29±0.05 


0.39±0.19 






5.00 


6.50 


± 


± 


0.12±0.03 


0.24±0.16 






6.50 


8.00 


± 


± 


± 


0.027±0.035 


0.150 


0.200 


0.50 


1.00 


0.18±0.06 


0.79±0.14 


1.35±0.19 


2.46±1.08 






1.00 


1.50 


0.26±0.06 


0.64±0.10 


0.83±0.11 


0.97±0.60 






1.50 


2.00 


0.31±0.06 


0.37±0.06 


0.61±0.09 








2.00 


2.50 


0.26±0.05 


0.40±0.07 


0.53±0.09 








2.50 


3.00 


± 


0.18±0.04 


0.30±0.07 


0.19±0.29 






3.00 


3.50 


± 


0.20±0.06 


0.26±0.05 


0.45±0.41 






3.50 


4.00 


± 


0.15±0.04 


0.12±0.03 


0.210±0.250 






4.00 


5.00 


± 


0.09±0.05 


0.08±0.02 


0.127±0.185 






5.00 


6.50 


± 


± 


0.031±0.013 


0.145±0.162 






6.50 


8.00 


± 


± 


± 


0.020±0.049 


0.200 


0.250 


0.50 


1.00 


0.26±0.14 


0.55±0.11 


0.86±0.15 


0.057±0.212 






1.00 


1.50 


0.11±0.04 


0.39±0.08 


0.42±0.09 


0.06±0.73 
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Tabic 31: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
cPa* /(dpdCl), for 7r+-Pb interactions at 3,5,8,12 GeV/c. Each row refers to a different (p m j n < p < 
Pmax, # m in — " < "max) bin, where p and 9 are the pion momentum and polar angle, respectively. The 
central value as well as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 32: HARP results for the double-differential tt + and tt~ production cross-section in the laboratory 
system, d 2 a* / (dpdfl) , for n~ - Be interactions at 8 and 12 GeV/c. Each row refers to a different (p m i n < 
P < Pmax, ^min < $ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. 
A finer angular binning than in the previous set of tables is used. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 
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Table 33: HARP results for the double-differential ir + and ir~ production cross-section in the laboratory 
system, d 2 a n / (dpdil) , for tt + - Be interactions at 8 and 8.9 GeV/c. Each row refers to a different 
(Pmin — P < Pmax, ^min — ^ < #max) bin, where p and are the pion momentum and polar angle, 
respectively. A finer angular binning than in the previous set of tables is used. The central value as well 
as the square-root of the diagonal elements of the covariancc matrix are given. 
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Table 34: HARP results for the double-differential ir + and tt~ production cross-section in the laboratory 
system, d 2 / (dpdQ) , for n~- C interactions at 8 and 12 GeV/c. Each row refers to a different (p m i n < 
P < Pmax, #min < $ < ^max) bin, where p and are the pion momentum and polar angle, respectively. 
A finer angular binning than in the previous set of tables is used. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 
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0.29±0.03 




2.50 


3.00 


0.10±0.02 


0.17±0.02 


0.22±0.02 


0.26±0.03 




3.00 


3.50 


0.09±0.02 


0.18±0.02 


0.22±0.02 


0.28±0.03 




3.50 


4.00 


0.090±0.013 


0.17±0.02 


0.20±0.02 


0.24±0.03 




4.00 


5.00 


0.083±0.010 


0.096±0.012 


0.21±0.02 


0.23±0.02 




5.00 


6.50 


0.045±0.006 


0.073±0.009 


0.17±0.02 


0.173±0.013 




6.50 


8.00 


± 


0.026±0.005 


± 


0.140±0.012 


0.075 0.100 


0.50 


0.75 


0.12±0.02 


0.14±0.03 


0.30±0.04 


0.37±0.06 




0.75 


1.00 


0.16±0.02 


0.23±0.03 


0.22±0.02 


0.22±0.03 




1.00 


1.25 


0.21±0.02 


0.27±0.03 


0.25±0.02 


0.29±0.03 




1.25 


1.50 


0.18±0.02 


0.30±0.03 


0.25±0.03 


0.33±0.04 




1.50 


2.00 


0.166±0.013 


0.25±0.02 


0.25±0.02 


0.30±0.03 




2.00 


2.50 


0.114±0.012 


0.21±0.02 


0.24±0.02 


0.34±0.03 




2.50 


3.00 


0.092±0.010 


0.18±0.02 


0.23±0.02 


0.32±0.03 




3.00 


3.50 


0.096±0.010 


0.17±0.02 


0.16±0.02 


0.25±0.02 




3.50 


4.00 


0.053±0.012 


0.136±0.014 


0.16±0.02 


0.23±0.02 




4.00 


5.00 


0.026±0.005 


0.072±0.008 


0.144±0.013 


0.171±0.014 




5.00 


6.50 


0.008±0.003 


0.039±0.006 


0.078±0.006 


0.105±0.008 




6.50 


8.00 


± 


0.008±0.002 


± 


0.052±0.006 
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Table 35: HARP results for the double-differential ir + and ir~ production cross-section in the laboratory 
system, d 2 er 71 " / (dpdfl) , for 7r+- C interactions at 8 GeV/c. Each row refers to a different {p m \ n < p < 
Pmax, ^min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of 
the diagonal elements of the covariancc matrix arc given. 



e ■ 


Q 

" lliclX 


^rnm 




d 2 a 7T+ /{dpdQ) 


d 2 a 7T /(dpdQ) 


(rad) 


(rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 










8 GeV/c 


8 GeV/c 


0.025 


0.050 


0.50 


0.75 


1 5+0 06 


27+0 08 






0.75 


1.00 


19+0 05 


19+0 05 






1.00 


1.25 


0.05±0.02 


0.12+0.03 






1.25 


1.50 


07±0 03 


18+0 06 

v. / . i v ' — i — v ' . v / w 






1.50 


2.00 


11±0 03 


17+0 04 






2.00 


2.50 


15±0 04 


12+0 03 






2.50 


3.00 


18+0 04 


09+0 02 






3.00 


3.50 


23+0 03 


15+0 05 






3.50 


4.00 


17+0 02 


10+0 02 






4.00 


5.00 


0.21+0.02 


0.11+0.02 






5.00 


6.50 


0.36+0.02 


0.097+0.014 


0.050 


0.075 


0.50 


0.75 


0.18+0.05 


0.16+0.04 






0.75 


1.00 


0.17+0.03 


0.12+0.03 






1.00 


1.25 


0.14+0.03 


0.17+0.04 






1.25 


1.50 


0.25+0.04 


0.15+0.03 






1.50 


2.00 


0.17+0.03 


0.15+0.02 






2.00 


2.50 


0.16+0.03 


0.16+0.02 






2.50 


3.00 


0.22+0.03 


0.11+0.02 






3.00 


3.50 


0.17+0.03 


0.081+0.014 






3.50 


4.00 


0.14+0.02 


0.13+0.02 






4.00 


5.00 


0.26+0.02 


0.065+0.010 






5.00 


6.50 


0.25+0.02 


0.046+0.006 


0.075 


0.100 


0.50 


0.75 


0.13+0.03 


0.17+0.04 






0.75 


1.00 


0.22+0.04 


0.13+0.02 






1.00 


1.25 


0.23+0.03 


0.14+0.03 






1.25 


1.50 


0.21+0.03 


0.15+0.02 






1.50 


2.00 


0.20+0.02 


0.15+0.02 






2.00 


2.50 


0.20+0.03 


0.11+0.02 






2.50 


3.00 


0.19+0.02 


0.12+0.02 






3.00 


3.50 


0.20+0.02 


0.081+0.011 






3.50 


4.00 


0.14+0.02 


0.088+0.013 






4.00 


5.00 


0.153+0.014 


0.039+0.005 






5.00 


6.50 


0.12+0.04 


0.011+0.002 
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Table 36: HARP results for the double-differential tt + and tt~ production cross-section in the laboratory 
system, d 2 er 71 " / (dpdCl) , for tt~- Al interactions at 8 and 12 GeV/c. Each row refers to a different (p m j n < 
P < Pmax, ^min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. 
A finer angular binning than in the previous set of tables is used. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 



^min $max 


Pmin 


Pmax 


d 2 a* + 


/{dpdQ) 




/(dpdQ) 


(rad) (rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 








8 GoV/c 


12 GeV/c 


8 GoV/c 


12 GeV/c 


0.025 0.050 


0.50 


0.75 


0.12±0.04 


0.32±0.09 


0.49±0.11 


0.50±0.12 




0.75 


1.00 


0.27±0.05 


0.37±0.06 


0.41±0.06 


0.42±0.07 




1.00 


1.25 


0.22±0.04 


0.33±0.06 


0.36±0.06 


0.38±0.06 




1.25 


1.50 


0.42±0.07 


0.38±0.06 


0.34±0.06 


0.50±0.09 




1.50 


2.00 


0.33±0.05 


0.36±0.06 


0.45±0.05 


0.51±0.06 




2.00 


2.50 


0.32±0.05 


0.48±0.06 


0.42±0.05 


0.49±0.07 




2.50 


3.00 


0.19±0.04 


0.45±0.06 


0.40±0.05 


0.37±0.05 




3.00 


3.50 


0.15±0.04 


0.30±0.04 


0.38±0.04 


0.43±0.05 




3.50 


4.00 


0.17±0.04 


0.31±0.04 


0.31±0.04 


0.41±0.04 




4.00 


5.00 


0.15±0.03 


0.26±0.03 


0.32±0.04 


0.41±0.04 




5.00 


6.50 


0.14±0.02 


0.19±0.02 


0.35±0.04 


0.39±0.03 




6.50 


8.00 


± 


0.16±0.02 


± 


0.40±0.03 


0.050 0.075 


0.50 


0.75 


0.28±0.06 


0.23±0.05 


0.55±0.09 


0.48±0.08 




0.75 


1.00 


0.25±0.04 


0.38±0.05 


0.45±0.05 


0.38±0.05 




1.00 


1.25 


0.34±0.04 


0.39±0.05 


0.32±0.04 


0.38±0.05 




1.25 


1.50 


0.30±0.04 


0.43±0.06 


0.33±0.05 


0.43±0.06 




1.50 


2.00 


0.33±0.03 


0.35±0.04 


0.46±0.04 


0.48±0.05 




2.00 


2.50 


0.20±0.03 


0.34±0.04 


0.40±0.04 


0.47±0.05 




2.50 


3.00 


0.13±0.03 


0.28±0.03 


0.39±0.04 


0.45±0.04 




3.00 


3.50 


0.18±0.03 


0.35±0.04 


0.33±0.04 


0.48±0.05 




3.50 


4.00 


0.17±0.02 


0.24±0.03 


0.33±0.04 


0.48±0.04 




4.00 


5.00 


0.12±0.02 


0.16±0.02 


0.33±0.03 


0.31±0.03 




5.00 


6.50 


0.073±0.010 


0.098±0.011 


0.29±0.02 


0.27±0.02 




6.50 


8.00 


± 


0.049±0.007 


± 


0.23±0.02 


0.075 0.100 


0.50 


0.75 


0.26±0.05 


0.25±0.04 


0.57±0.08 


0.63±0.09 




0.75 


1.00 


0.35±0.04 


0.40±0.05 


0.43±0.04 


0.43±0.05 




1.00 


1.25 


0.36±0.04 


0.52±0.05 


0.44±0.04 


0.49±0.05 




1.25 


1.50 


0.40±0.04 


0.42±0.04 


0.50±0.05 


0.52±0.06 




1.50 


2.00 


0.28±0.02 


0.39±0.03 


0.41±0.03 


0.54±0.04 




2.00 


2.50 


0.23±0.02 


0.31±0.03 


0.44±0.04 


0.48±0.04 




2.50 


3.00 


0.14±0.02 


0.30±0.02 


0.40±0.03 


0.48±0.04 




3.00 


3.50 


0.17±0.02 


0.25±0.02 


0.28±0.03 


0.47±0.03 




3.50 


4.00 


0.10±0.02 


0.21±0.02 


0.24±0.03 


0.38±0.03 




4.00 


5.00 


0.060±0.009 


0.124±0.012 


0.23±0.02 


0.23±0.02 




5.00 


6.50 


0.015±0.005 


0.052±0.006 


0.140±0.011 


0.146±0.011 




6.50 


8.00 


± 


0.016±0.003 


± 


0.071±0.007 
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Table 37: HARP results for the double-differential tt + and tt~ production cross-section in the laboratory 
system, d 2 cr*" / (dpdft) , for 7r + - Al interactions at 8 and 12.9 GeV/c. Each row refers to a different 
(Pmin < P < Pmax, ^min — ^ < ^max) bin, where p and 9 are the pion momentum and polar angle, 
respectively. A finer angular binning than in the previous set of tables is used. The central value as well 
as the square-root of the diagonal elements of the covariancc matrix are given. 



^min 


#max 


Pmin 


Pmax 


d 2 a" + 


/(dpdQ.) 


d 2 0^ 


/{dpdQ) 


(rad) 


(rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 










8 GcV/c 


12.9 GeV/c 


8 GcV/c 


12.9 GoV/c 


0.025 


0.050 


0.50 


0.75 


0.37±0.11 


0.51±0.16 


0.15±0.06 


0.35±0.16 






0.75 


1.00 


0.38±0.08 


0.43±0.13 


0.37±0.09 


0.38±0.12 






1.00 


1.25 


0.24±0.07 


0.36±0.12 


0.27±0.07 


0.36±0.14 






1.25 


1.50 


0.11±0.04 


0.27±0.10 


0.48±0.12 


0.13±0.08 






1.50 


2.00 


0.23±0.05 


0.54±0.10 


0.26±0.06 


0.30±0.09 






2.00 


2.50 


0.27±0.06 


0.42±0.09 


0.24±0.05 


0.34±0.11 






2.50 


3.00 


0.35±0.06 


0.63±0.10 


0.21±0.04 


0.27±0.07 






3.00 


3.50 


0.31±0.05 


0.47±0.08 


0.25±0.04 


0.35±0.08 






3.50 


4.00 


0.38±0.05 


0.49±0.09 


0.21±0.04 


0.40±0.08 






4.00 


5.00 


0.37±0.04 


0.39±0.06 


0.19±0.03 


0.21±0.04 






5.00 


6.50 


0.50±0.03 


0.41±0.05 


0.19±0.03 


0.14±0.03 






6.50 


8.00 


± 


0.42±0.05 


± 


0.16±0.03 


0.050 


0.075 


0.50 


0.75 


0.38±0.09 


0.48±0.11 


0.22±0.06 


0.21±0.08 






0.75 


1.00 


0.33±0.06 


0.64±0.10 


0.15±0.04 


0.27±0.07 






1.00 


1.25 


0.33±0.06 


0.36±0.08 


0.24±0.05 


0.40±0.09 






1.25 


1.50 


0.21±0.05 


0.36±0.08 


0.18±0.04 


0.39±0.09 






1.50 


2.00 


0.28±0.04 


0.43±0.07 


0.23±0.04 


0.43±0.06 






2.00 


2.50 


0.37±0.05 


0.52±0.07 


0.24±0.04 


0.33±0.06 






2.50 


3.00 


0.34±0.05 


0.45±0.06 


0.16±0.03 


0.24±0.05 






3.00 


3.50 


0.33±0.05 


0.46±0.06 


0.21±0.03 


0.32±0.06 






3.50 


4.00 


0.29±0.04 


0.42±0.06 


0.18±0.03 


0.29±0.05 






4.00 


5.00 


0.42±0.04 


0.41±0.04 


0.13±0.02 


0.17±0.03 






5.00 


6.50 


0.40±0.03 


0.29±0.03 


0.080±0.011 


0.14±0.02 






6.50 


8.00 


± 


0.21±0.03 


± 


0.052±0.012 


0.075 


0.100 


0.50 


0.75 


0.26±0.06 


0.35±0.09 


0.32±0.07 


0.24±0.08 






0.75 


1.00 


0.26±0.05 


0.59±0.08 


0.26±0.04 


0.30±0.07 






1.00 


1.25 


0.34±0.05 


0.47±0.07 


0.33±0.05 


0.41±0.07 






1.25 


1.50 


0.40±0.06 


0.45±0.07 


0.28±0.04 


0.41±0.07 






1.50 


2.00 


0.35±0.04 


0.43±0.06 


0.26±0.03 


0.36±0.05 






2.00 


2.50 


0.34±0.04 


0.54±0.06 


0.27±0.03 


0.36±0.04 






2.50 


3.00 


0.33±0.04 


0.34±0.04 


0.15±0.02 


0.29±0.04 






3.00 


3.50 


0.36±0.04 


0.48±0.05 


0.12±0.02 


0.27±0.03 






3.50 


4.00 


0.31±0.05 


0.36±0.04 


0.09±0.02 


0.29±0.03 






4.00 


5.00 


0.31±0.02 


0.25±0.03 


0.043±0.007 


0.12±0.02 






5.00 


6.50 


0.16±0.05 


0.14±0.02 


0.025±0.005 


0.051±0.010 






6.50 


8.00 


± 


0.080±0.012 


± 


0.015±0.004 
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Table 38: HARP results for the double-differential tt + and tt~ production cross-section in the laboratory 
system, cPa* / (dpdCt), for tt~- Cu interactions at 8 and 12 GeV/c. Each row refers to a different (p m i n < 
P < Pmax, ^min < $ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. 
A finer angular binning than in the previous set of tables is used. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 



^min $max 


Pmin 


Pmax 


d 2 a* + 


/{dpdQ) 


d 2 ^' 


/(dpdQ) 


(rad) (rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 








8 GoV/c 


12 GeV/c 


8 GcV/c 


12 GeV/c 


0.025 0.050 


0.50 


0.75 


0.44±0.10 


0.73±0.19 


1.03±0.16 


1.23±0.27 




0.75 


1.00 


0.47±0.06 


0.68±0.12 


0.55±0.07 


1.08±0.16 




1.00 


1.25 


0.42±0.05 


0.64±0.11 


0.46±0.06 


0.54±0.10 




1.25 


1.50 


0.70±0.09 


0.80±0.14 


0.55±0.09 


0.89±0.17 




1.50 


2.00 


0.55±0.07 


0.71±0.11 


0.78±0.07 


0.93±0.12 




2.00 


2.50 


0.50±0.07 


0.54±0.09 


0.68±0.06 


0.80±0.11 




2.50 


3.00 


0.46±0.07 


0.73±0.11 


0.57±0.06 


0.73±0.10 




3.00 


3.50 


0.24±0.06 


0.76±0.11 


0.42±0.05 


0.79±0.10 




3.50 


4.00 


0.24±0.05 


0.61±0.09 


0.35±0.05 


0.63±0.08 




4.00 


5.00 


0.29±0.05 


0.33±0.05 


0.37±0.05 


0.61±0.07 




5.00 


6.50 


0.24±0.04 


0.34±0.04 


0.41±0.05 


0.57±0.06 




6.50 


8.00 


± 


0.27±0.04 


± 


0.61±0.05 


0.050 0.075 


0.50 


0.75 


0.55±0.09 


0.81±0.15 


1.12±0.15 


1.25±0.22 




0.75 


1.00 


0.63±0.07 


0.72±0.11 


0.77±0.07 


1.09±0.14 




1.00 


1.25 


0.59±0.06 


0.85±0.12 


0.61±0.06 


0.99±0.13 




1.25 


1.50 


0.67±0.07 


0.83±0.12 


0.62±0.07 


0.81±0.11 




1.50 


2.00 


0.61±0.05 


0.63±0.08 


0.76±0.06 


0.92±0.10 




2.00 


2.50 


0.37±0.05 


0.60±0.07 


0.70±0.06 


0.93±0.10 




2.50 


3.00 


0.24±0.04 


0.49±0.06 


0.62±0.05 


0.67±0.08 




3.00 


3.50 


0.29±0.05 


0.56±0.07 


0.50±0.06 


0.70±0.08 




3.50 


4.00 


0.25±0.04 


0.47±0.07 


0.43±0.06 


0.63±0.09 




4.00 


5.00 


0.23±0.03 


0.28±0.04 


0.55±0.06 


0.59±0.05 




5.00 


6.50 


0.12±0.02 


0.17±0.02 


0.43±0.04 


0.48±0.04 




6.50 


8.00 


± 


0.12±0.02 


± 


0.34±0.03 


0.075 0.100 


0.50 


0.75 


0.55±0.09 


0.72±0.12 


1.08±0.13 


1.42±0.19 




0.75 


1.00 


0.75±0.07 


0.59±0.08 


0.70±0.06 


0.74±0.09 




1.00 


1.25 


0.67±0.06 


0.87±0.10 


0.82±0.07 


1.10±0.12 




1.25 


1.50 


0.58±0.05 


0.65±0.07 


0.75±0.07 


0.89±0.10 




1.50 


2.00 


0.48±0.04 


0.73±0.08 


0.74±0.05 


0.96±0.08 




2.00 


2.50 


0.40±0.04 


0.75±0.07 


0.67±0.05 


1.10±0.09 




2.50 


3.00 


0.29±0.03 


0.58±0.05 


0.63±0.05 


0.83±0.07 




3.00 


3.50 


0.30±0.03 


0.51±0.05 


0.45±0.05 


0.84±0.08 




3.50 


4.00 


0.16±0.04 


0.35±0.04 


0.39±0.04 


0.71±0.08 




4.00 


5.00 


0.082±0.014 


0.24±0.03 


0.33±0.03 


0.40±0.03 




5.00 


6.50 


0.024±0.007 


0.089±0.014 


0.21±0.02 


0.24±0.02 




6.50 


8.00 


± 


0.018±0.005 


± 


0.122±0.013 
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Table 39: HARP results for the double-differential ir + and ir~ production cross-section in the laboratory 
system, d 2 er 7r / (dpdft) , for 7r + - Cu interactions at 8 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, ^min — ^ < ^max) bin, where p and are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of 
the diagonal elements of the covariancc matrix arc given. 



e ■ 


Q 

" lliclX 






d 2 a 7T+ /{dpdQ) 


d 2 a 7T /(dpdQ) 


(rad) 


(rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 










8 GeV/c 


8 GeV/c 


0.025 


0.050 


0.50 


0.75 


22+0 08 


45+0 1 5 






0.75 


1.00 


32±0 08 


45+0 11 

V ' . 1 '7 — 1 — V / . 1 1 






1.00 


1.25 


0.26±0.07 


0.37+0.10 






1.25 


1.50 


0.49±0.12 


0.31+0.11 






1.50 


2.00 


32±0 08 


0.64+0.12 






2.00 


2.50 


34±0 09 


36+0 08 






2.50 


3.00 


61 +0 1 2 


27+0 06 






3.00 


3.50 


32±0 06 


40+0 07 






3.50 


4.00 


47±0 07 


31+0 05 






4.00 


5.00 


0.40±0.05 


0.32+0.04 






5.00 


6.50 


0.76±0.05 


0.30+0.04 


0.050 


0.075 


0.50 


0.75 


0.57±0.14 


0.44+0.11 






0.75 


1.00 


0.59±0.11 


0.49+0.09 






1.00 


1.25 


0.49±0.09 


0.48+0.10 






1.25 


1.50 


0.66±0.11 


0.56+0.10 






1.50 


2.00 


0.60±0.08 


0.44+0.07 






2.00 


2.50 


0.59±0.09 


0.36+0.06 






2.50 


3.00 


0.51±0.08 


0.31+0.05 






3.00 


3.50 


0.45±0.08 


0.30+0.05 






3.50 


4.00 


0.44±0.07 


0.22+0.04 






4.00 


5.00 


0.58±0.06 


0.23+0.03 






5.00 


6.50 


0.62±0.04 


0.12+0.02 


0.075 


0.100 


0.50 


0.75 


0.64±0.15 


0.61+0.12 






0.75 


1.00 


0.53±0.09 


0.47+0.08 






1.00 


1.25 


0.66±0.10 


0.46+0.08 






1.25 


1.50 


0.59±0.08 


0.64+0.09 






1.50 


2.00 


0.55±0.07 


0.39+0.05 






2.00 


2.50 


0.54±0.07 


0.40+0.05 






2.50 


3.00 


0.56±0.07 


0.24+0.03 






3.00 


3.50 


0.41±0.05 


0.22+0.03 






3.50 


4.00 


0.47±0.07 


0.15+0.03 






4.00 


5.00 


0.41±0.04 


0.13+0.02 






5.00 


6.50 


0.27±0.09 


0.024+0.005 
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Table 40: HARP results for the double-differential ir + and 7r~ production cross-section in the laboratory 
system, d 2 a* / (dpdCt) , for tt~- Sn interactions at 8 and 12 GeV/c. Each row refers to a different (p m i n < 
P < Pmax, #min < $ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. 
A finer angular binning than in the previous set of tables is used. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 



^min $max 


Pmin 


Pmax 


d 2 a* + 


/{dpdQ) 


d 2 a^ 


/(dpdQ,) 


(rad) (rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 








8 GoV/c 


12 GeV/c 


8 GeV/c 


12 GeV/c 


0.025 0.050 


0.50 


0.75 


0.64±0.14 


1.56±0.30 


1.39±0.21 


1.86±0.34 




0.75 


1.00 


0.53±0.07 


1.16±0.19 


0.70±0.08 


1.53±0.17 




1.00 


1.25 


0.36±0.05 


0.94±0.13 


0.29±0.04 


0.93±0.12 




1.25 


1.50 


0.55±0.09 


1.15±0.17 


0.55±0.10 


1.34±0.20 




1.50 


2.00 


0.63±0.10 


1.01±0.76 


0.80±0.08 


1.38±0.15 




2.00 


2.50 


0.84±0.12 


1.07±0.36 


0.75±0.09 


1.19±0.16 




2.50 


3.00 


0.53±0.11 


0.88±0.15 


0.70±0.09 


1.08±0.12 




3.00 


3.50 


0.24±0.08 


0.92±0.16 


0.45±0.07 


1.04±0.11 




3.50 


4.00 


0.37±0.09 


0.81±0.10 


0.27±0.05 


0.95±0.13 




4.00 


5.00 


0.36±0.07 


0.58±0.07 


0.43±0.07 


0.91±0.09 




5.00 


6.50 


0.25±0.05 


0.43±0.05 


0.63±0.09 


0.73±0.07 




6.50 


8.00 


± 


0.39±0.05 


± 


0.84±0.07 


0.050 0.075 


0.50 


0.75 


0.55±0.11 


0.84±0.16 


1.19±0.17 


1.98±0.28 




0.75 


1.00 


0.57±0.09 


0.99±0.13 


l.OliO.ll 


1.13±0.14 




1.00 


1.25 


0.79±0.10 


1.01±0.13 


0.82±0.10 


1.20±0.14 




1.25 


1.50 


0.96±0.12 


1.24±0.17 


0.88±0.12 


1.17±0.14 




1.50 


2.00 


0.77±0.07 


l.OliO.ll 


0.93±0.09 


1.50±0.14 




2.00 


2.50 


0.58±0.08 


0.95±0.11 


0.84±0.09 


1.20±0.11 




2.50 


3.00 


0.32±0.07 


0.68±0.07 


0.69±0.08 


1.19±0.11 




3.00 


3.50 


0.39±0.08 


0.87±0.10 


0.82±0.11 


1.10±0.11 




3.50 


4.00 


0.31±0.06 


0.60±0.07 


0.58±0.09 


1.03±0.11 




4.00 


5.00 


0.23±0.04 


0.41±0.05 


0.70±0.09 


0.79±0.07 




5.00 


6.50 


0.16±0.03 


0.26±0.03 


0.56±0.06 


0.67±0.05 




6.50 


8.00 


± 


0.14±0.02 


± 


0.45±0.04 


0.075 0.100 


0.50 


0.75 


0.73±0.12 


1.09±0.16 


1.61±0.20 


2.13±0.26 




0.75 


1.00 


0.80±0.10 


1.54±0.16 


1.13±0.11 


1.56±0.15 




1.00 


1.25 


1.08±0.11 


1.37±0.13 


0.98±0.09 


1.31±0.12 




1.25 


1.50 


0.73±0.08 


1.30±0.11 


0.88±0.10 


1.43±0.15 




1.50 


2.00 


0.68±0.06 


1.03±0.09 


0.94±0.08 


1.29±0.10 




2.00 


2.50 


0.48±0.05 


0.92±0.07 


0.96±0.09 


1.43±0.11 




2.50 


3.00 


0.41±0.05 


0.80±0.07 


0.93±0.08 


1.28±0.10 




3.00 


3.50 


0.34±0.04 


0.56±0.05 


0.54±0.08 


0.97±0.07 




3.50 


4.00 


0.14±0.05 


0.43±0.05 


0.51±0.07 


0.88±0.07 




4.00 


5.00 


0.09±0.02 


0.33±0.03 


0.47±0.05 


0.57±0.06 




5.00 


6.50 


0.021±0.008 


0.14±0.02 


0.26±0.03 


0.31±0.02 




6.50 


8.00 


± 


0.037±0.008 


± 


0.17±0.02 



74 



Table 41: HARP results for the double-differential ir + and ir~ production cross-section in the laboratory 
system, d 2 a v / (dpdfl) , for n + - Sn interactions at 8 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, # m in — ^ < "max) bin, where p and are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of 
the diagonal elements of the covariancc matrix arc given. 



e ■ 


a 

" IllcLX 






d 2 a 7T+ /{dpdQ) 


d 2 a 7T /(dpdQ) 


(rad) 


(rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 










8 GeV/c 


8 GeV/c 


0.025 


0.050 


0.50 


0.75 


78+0 18 


58+0 1 4 






0.75 


1.00 


63±0 10 


57±0 11 






1.00 


1.25 


0.34±0.07 


0.68+0.14 






1.25 


1.50 


27±0 08 


0.34+0.11 






1.50 


2.00 


63±0 13 


82+0 16 






2.00 


2.50 


55±0 13 


59+0 11 






2.50 


3.00 


5Q+0 1 3 


67+0 1 1 






3.00 


3.50 


52±0 09 


48+0 08 






3.50 


4.00 


54±0 07 


45+0 07 






4.00 


5.00 


0.61±0.06 


0.52+0.07 






5.00 


6.50 


0.98±0.06 


0.39+0.06 


0.050 


0.075 


0.50 


0.75 


0.85±0.19 


0.88+0.19 






0.75 


1.00 


0.61±0.12 


0.54+0.11 






1.00 


1.25 


0.60±0.11 


0.66+0.13 






1.25 


1.50 


0.97±0.16 


0.81+0.15 






1.50 


2.00 


0.85±0.11 


0.49+0.08 






2.00 


2.50 


0.67±0.11 


0.61+0.10 






2.50 


3.00 


0.54±0.09 


0.41+0.07 






3.00 


3.50 


0.44±0.09 


0.42+0.07 






3.50 


4.00 


0.39±0.08 


0.39+0.06 






4.00 


5.00 


0.77±0.08 


0.28+0.04 






5.00 


6.50 


0.88±0.06 


0.14+0.02 


0.075 


0.100 


0.50 


0.75 


1.33±0.26 


0.91+0.17 






0.75 


1.00 


1.01±0.15 


0.74+0.11 






1.00 


1.25 


0.82±0.12 


0.67+0.11 






1.25 


1.50 


0.79±0.12 


0.68+0.10 






1.50 


2.00 


0.70±0.23 


0.69+0.08 






2.00 


2.50 


0.68±0.35 


0.44+0.06 






2.50 


3.00 


0.72±0.28 


0.40+0.05 






3.00 


3.50 


0.62±0.08 


0.32+0.05 






3.50 


4.00 


0.65±0.11 


0.22+0.03 






4.00 


5.00 


0.62±0.08 


0.13+0.02 






5.00 


6.50 


0.42±0.15 


0.044+0.009 
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Table 42: HARP results for the double-differential ir + and tt~ production cross-section in the laboratory 
system, d 2 a* / (dpdCt) , for tt~- Ta interactions at 8 and 12 GeV/c. Each row refers to a different (p m i n < 
P < Pmax, #min < $ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. 
A finer angular binning than in the previous set of tables is used. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 



^min $max 


Pmin 


Pmax 


d 2 a* + 


/{dpdQ) 


d 2 a^ 


/(dpdQ,) 


(rad) (rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 








8 GoV/c 


12 GeV/c 


8 GeV/c 


12 GeV/c 


0.025 0.050 


0.50 


0.75 


1.47±0.25 


1.65±0.32 


2.65±0.39 


2.44±0.39 




0.75 


1.00 


0.72±0.08 


1.11±0.15 


1.18±0.12 


1.34±0.16 




1.00 


1.25 


0.57±0.07 


1.02±0.15 


0.43±0.06 


0.89±0.13 




1.25 


1.50 


0.90±0.12 


1.19±0.18 


0.76±0.13 


1.61±0.24 




1.50 


2.00 


0.72±0.11 


1.21±0.17 


1.09±0.11 


1.51±0.19 




2.00 


2.50 


0.90±0.15 


1.26±0.17 


1.10±0.12 


1.80±0.23 




2.50 


3.00 


0.90±0.15 


1.26±0.16 


0.92±0.13 


1.34±0.18 




3.00 


3.50 


0.51±0.12 


1.49±0.18 


0.50±0.08 


1.50±0.16 




3.50 


4.00 


0.51±0.11 


0.77±0.14 


0.32±0.06 


1.08±0.13 




4.00 


5.00 


0.44±0.08 


0.55±0.08 


0.44±0.08 


0.93±0.10 




5.00 


6.50 


0.35±0.07 


0.55±0.07 


0.53±0.09 


0.92±0.11 




6.50 


8.00 


± 


0.36±0.05 


± 


0.96±0.09 


0.050 0.075 


0.50 


0.75 


0.82±0.14 


1.39±0.25 


1.27±0.18 


2.13±0.31 




0.75 


1.00 


0.75±0.11 


1.38±0.18 


1.12±0.12 


1.35±0.18 




1.00 


1.25 


1.12±0.14 


1.27±0.17 


0.89±0.11 


1.62±0.20 




1.25 


1.50 


1.41±0.16 


1.87±0.25 


1.08±0.15 


1.92±0.22 




1.50 


2.00 


1.21±0.11 


1.23±0.13 


1.43±0.12 


1.64±0.16 




2.00 


2.50 


0.77±0.11 


1.16±0.12 


1.21±0.13 


1.70±0.15 




2.50 


3.00 


0.32±0.09 


1.02±0.11 


l.lliO.ll 


1.43±0.14 




3.00 


3.50 


0.40±0.10 


0.95±0.11 


0.89±0.13 


1.41±0.13 




3.50 


4.00 


0.42±0.08 


0.87±0.12 


0.82±0.13 


1.06±0.11 




4.00 


5.00 


0.24±0.04 


0.54±0.06 


0.75±0.10 


0.88±0.08 




5.00 


6.50 


0.21±0.04 


0.31±0.04 


0.73±0.09 


0.73±0.06 




6.50 


8.00 


± 


0.19±0.03 


± 


0.53±0.05 


0.075 0.100 


0.50 


0.75 


1.21±0.19 


1.42±0.22 


2.19±0.26 


2.67±0.34 




0.75 


1.00 


1.31±0.14 


2.10±0.22 


1.17±0.12 


2.23±0.21 




1.00 


1.25 


1.27±0.12 


1.77±0.17 


1.19±0.12 


1.78±0.17 




1.25 


1.50 


1.27±0.12 


1.71±0.16 


1.16±0.13 


1.86±0.20 




1.50 


2.00 


0.93±0.08 


1.30±0.11 


1.06±0.09 


1.62±0.13 




2.00 


2.50 


0.66±0.07 


1.03±0.10 


1.24±0.11 


1.51±0.13 




2.50 


3.00 


0.39±0.06 


0.79±0.08 


0.92±0.08 


1.53±0.12 




3.00 


3.50 


0.35±0.05 


0.84±0.08 


0.75±0.09 


1.17±0.09 




3.50 


4.00 


0.19±0.06 


0.57±0.06 


0.61±0.07 


0.98±0.08 




4.00 


5.00 


0.14±0.03 


0.35±0.04 


0.61±0.07 


0.69±0.06 




5.00 


6.50 


0.032±0.013 


0.15±0.02 


0.31±0.04 


0.36±0.03 




6.50 


8.00 


± 


0.040±0.009 


± 


0.17±0.02 
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Table 43: HARP results for the double-differential ir + and ir~ production cross-section in the laboratory 
system, d 2 cr 71 " / '(dpdfl) , for 7r + - Ta interactions at 8 GeV/c. Each row refers to a different (p m i n < p < 
Pmax, 6* m i n <0< #max) bin, where p and 9 are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of 
the diagonal elements of the covariancc matrix arc given. 



e ■ 


a 

" ILLcLX 






d 2 a 7T+ /(dpdQ.) 


d 2 a^/{dpdQ.) 


(rad) 


(rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 










8 GeV/c 


8 GeV/c 


0.025 


0.050 


0.50 


0.75 


1 66+0 37 


92+0 25 






0.75 


1.00 


1 02±0 16 

1 . V ' — — 1 — VI. 1 \ 1 


60+0 13 






1.00 


1.25 


0.50±0.10 


0.96+0.21 






1.25 


1.50 


39+0 12 

V / . ' 1 - J — 1 — Uili 


33+0 13 






1.50 


2.00 


66±0 16 


87+0 19 






2.00 


2.50 


62±0 17 


63+0 14 






2.50 


3.00 


72+0 18 


54+0 1 1 






3.00 


3.50 


58±0 12 

\J m 'J \J 1 Utli 


82+0 15 

\J m yj £j 1 \J m _L \J 






3.50 


4.00 


59±0 11 


63+0 11 

V / . vl ' 1 — 1 — V 1 . 1 1 






4.00 


5.00 


0.56±0.08 


0.56+0.09 






5.00 


6.50 


1.12±0.08 


0.48+0.07 


0.050 


0.075 


0.50 


0.75 


1.15±0.27 


0.79+0.22 






0.75 


1.00 


0.90±0.19 


0.72+0.16 






1.00 


1.25 


0.95±0.19 


1.02+0.22 






1.25 


1.50 


0.96±0.20 


0.54+0.13 






1.50 


2.00 


0.57±0.12 


0.98+0.16 






2.00 


2.50 


0.70±0.15 


0.60+0.11 






2.50 


3.00 


0.72±0.14 


0.62+0.12 






3.00 


3.50 


0.46±0.12 


0.36+0.07 






3.50 


4.00 


0.42±0.10 


0.41+0.08 






4.00 


5.00 


0.95±0.11 


0.31+0.05 






5.00 


6.50 


1.06±0.08 


0.18+0.03 


0.075 


0.100 


0.50 


0.75 


1.11±0.27 


0.88+0.21 






0.75 


1.00 


1.04±0.18 


0.75+0.14 






1.00 


1.25 


1.15±0.19 


1.19+0.20 






1.25 


1.50 


1.03±0.17 


0.79+0.14 






1.50 


2.00 


0.67±0.11 


0.62+0.10 






2.00 


2.50 


0.72±0.12 


0.47+0.08 






2.50 


3.00 


0.74±0.12 


0.51+0.08 






3.00 


3.50 


0.72±0.10 


0.44+0.07 






3.50 


4.00 


0.70±0.11 


0.22+0.04 






4.00 


5.00 


0.71±0.07 


0.19+0.03 






5.00 


6.50 


0.41±0.15 


0.051+0.011 
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Table 44: HARP results for the double-differential ir + and tt~ production cross-section in the laboratory 
system, cPa* / (dpdQ), for tt~- Pb interactions at 8 and 12 GeV/c. Each row refers to a different (p m i n < 
P < Pmax, ^min < $ < ^max) bin, where p and 9 are the pion momentum and polar angle, respectively. 
A finer angular binning than in the previous set of tables is used. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 



^min $max 


Pmin 


Pmax 


d 2 a* + 


/(dpdQ) 


d 2 a^ 


/(dpdQ,) 


(rad) (rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GcV/c rad)) 








8 GoV/c 


12 GeV/c 


8 GeV/c 


12 GeV/c 


0.025 0.050 


0.50 


0.75 


1.60±0.27 


1.58±0.29 


3.24±0.43 


2.52±0.41 




0.75 


1.00 


0.87±0.09 


1.37±0.16 


1.26±0.14 


1.74±0.18 




1.00 


1.25 


0.56±0.07 


1.06±0.14 


0.54±0.07 


0.88±0.13 




1.25 


1.50 


0.76±0.11 


1.25±0.18 


0.54±0.09 


0.84±0.14 




1.50 


2.00 


0.82±0.13 


1.40±0.18 


1.27±0.14 


1.53±0.20 




2.00 


2.50 


1.12±0.16 


1.24±0.17 


1.24±0.14 


1.66±0.20 




2.50 


3.00 


0.76±0.15 


1.28±0.16 


1.02±0.13 


0.87±0.13 




3.00 


3.50 


0.41±0.12 


0.92±0.14 


0.63±0.10 


1.27±0.16 




3.50 


4.00 


0.53±0.12 


0.85±0.13 


0.42±0.09 


1.05±0.13 




4.00 


5.00 


0.47±0.09 


0.56±0.09 


0.54±0.10 


0.73±0.10 




5.00 


6.50 


0.40±0.08 


0.44±0.07 


0.65±0.11 


0.62±0.09 




6.50 


8.00 


± 


0.33±0.06 


± 


0.77±0.10 


0.050 0.075 


0.50 


0.75 


0.92±0.15 


1.45±0.22 


1.91±0.24 


2.46±0.33 




0.75 


1.00 


0.75±0.11 


1.33±0.17 


1.39±0.15 


1.34±0.18 




1.00 


1.25 


1.00±0.13 


1.82±0.20 


1.12±0.14 


1.56±0.19 




1.25 


1.50 


1.72±0.20 


1.70±0.20 


1.15±0.16 


1.49±0.18 




1.50 


2.00 


1.19±0.10 


1.38±0.14 


1.43±0.14 


1.66±0.16 




2.00 


2.50 


0.71±0.10 


1.21±0.12 


1.24±0.13 


1.52±0.15 




2.50 


3.00 


0.36±0.10 


l.OOiO.10 


0.95±0.11 


1.29±0.13 




3.00 


3.50 


0.41±0.11 


1.04±0.11 


0.90±0.14 


1.16±0.12 




3.50 


4.00 


0.42±0.08 


0.76±0.10 


0.84±0.13 


0.93±0.11 




4.00 


5.00 


0.36±0.06 


0.43±0.06 


0.82±0.12 


0.76±0.08 




5.00 


6.50 


0.20±0.04 


0.29±0.04 


0.67±0.08 


0.65±0.06 




6.50 


8.00 


± 


0.13±0.02 


± 


0.51±0.05 


0.075 0.100 


0.50 


0.75 


0.96±0.16 


1.81±0.23 


2.21±0.26 


2.93±0.34 




0.75 


1.00 


1.25±0.15 


1.83±0.18 


1.28±0.13 


1.92±0.19 




1.00 


1.25 


1.75±0.16 


1.77±0.16 


1.27±0.12 


1.59±0.15 




1.25 


1.50 


1.30±0.13 


1.71±0.14 


1.53±0.18 


1.71±0.19 




1.50 


2.00 


1.00±0.09 


1.42±0.12 


1.38±0.11 


1.73±0.13 




2.00 


2.50 


0.61±0.07 


1.17±0.10 


1.14±0.10 


1.67±0.12 




2.50 


3.00 


0.45±0.06 


0.90±0.08 


1.04±0.10 


1.38±0.11 




3.00 


3.50 


0.43±0.06 


0.91±0.08 


0.74±0.11 


1.22±0.09 




3.50 


4.00 


0.24±0.06 


0.55±0.07 


0.63±0.09 


0.99±0.10 




4.00 


5.00 


0.12±0.03 


0.32±0.04 


0.46±0.07 


0.54±0.06 




5.00 


6.50 


0.033±0.014 


0.14±0.02 


0.32±0.04 


0.33±0.03 




6.50 


8.00 


± 


0.032±0.008 


± 


0.19±0.02 



78 



Table 45: HARP results for the double-differential ir + and ir~ production cross-section in the laboratory 
system, d 2 er 7r / (dpdCt) , for tt + - Pb interactions at 8 GeV/c. Each row refers to a different (p m { n < p < 
Pmax, ^min — ^ < "max) bin, where p and are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of 
the diagonal elements of the covariancc matrix arc given. 



e ■ 


a 

" ILLcLX 






d 2 a 7T+ /{dpdQ) 


d 2 a 7T /(dpdQ) 


(rad) 


(rad) 


(GeV/c) 


(GeV/c) 


(barn/(GeV/c rad)) 


(barn/(GeV/c rad)) 










8 GeV/c 


8 GeV/c 


0.025 


0.050 


0.50 


0.75 


1 72+0 38 


75+0 21 






0.75 


1.00 


1 04±0 16 


85+0 17 






1.00 


1.25 


0.42±0.09 


0.73+0.17 






1.25 


1.50 


58±0 15 


79+0 25 






1.50 


2.00 


73±0 17 


66+0 17 






2.00 


2.50 


63±0 20 


72+0 16 






2.50 


3.00 


4Q+0 1 5 


60+0 1 3 






3.00 


3.50 


65±0 13 


51+0 10 






3.50 


4.00 


59±0 10 

V ' . '1-1 — 1 — V ' . 1 \ f 


67+0 12 

\ / . \> I — i — yj • _i_ £j 






4.00 


5.00 


0.63±0.08 


0.58+0.09 






5.00 


6.50 


1.18±0.08 


0.49+0.08 


0.050 


0.075 


0.50 


0.75 


0.66±0.18 


0.75+0.20 






0.75 


1.00 


0.77±0.18 


0.73+0.16 






1.00 


1.25 


0.43±0.12 


0.69+0.17 






1.25 


1.50 


0.75±0.18 


1.00+0.20 






1.50 


2.00 


0.69±0.13 


0.77+0.13 






2.00 


2.50 


0.74±0.15 


0.69+0.13 






2.50 


3.00 


0.85±0.16 


0.49+0.10 






3.00 


3.50 


0.55±0.13 


0.51+0.10 






3.50 


4.00 


0.50±0.11 


0.52+0.09 






4.00 


5.00 


1.09±0.12 


0.40+0.06 






5.00 


6.50 


1.11±0.09 


0.22+0.03 


0.075 


0.100 


0.50 


0.75 


1.04±0.27 


1.05+0.23 






0.75 


1.00 


1.02±0.19 


1.27+0.20 






1.00 


1.25 


0.86±0.16 


1.30+0.21 






1.25 


1.50 


0.99±0.17 


0.91+0.15 






1.50 


2.00 


0.84±0.13 


0.81+0.12 






2.00 


2.50 


0.88±0.15 


0.33+0.06 






2.50 


3.00 


0.66±0.11 


0.53+0.08 






3.00 


3.50 


0.66±0.09 


0.34+0.06 






3.50 


4.00 


0.82±0.13 


0.30+0.06 






4.00 


5.00 


0.75±0.07 


0.21+0.03 






5.00 


6.50 


0.44±0.15 


0.065+0.013 
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Normalization factors, incident tt 



+ 



o 
o 

(N 

X> 

<D 

PL. 
(N 



X 

in: 

CD 



> 

in 
o 

(N 

(N 
O 
On 
O 



X 



Generator 


Be 3 GeV 


Ta 3 GeV 


Be 5 GeV 


Ta 5 GeV 


Be 8 GeV 


Ta 8 GeV 


Be 1 2 GeV 




7T+ 


7T 


7T+ 


7T 


7T+ 


7T 


7T+ 


7T 


7T+ 


7T 


7T+ 


7T 


7T+ 


7T 


bert 


1.1 


1.1 


0.9 


0.6 


1.3 


1.3 


1.1 


1.1 


1.7 


1.4 


1.5 


1.4 


2.0 


1.4 


Ihep 


1.3 


0.6 


2.8 


1.0 


1.5 


0.8 


2.5 


1.1 


1.3 


0.8 


1.0 


0.5 


1.0 


0.8 


qgsb 


4.8 


1.8 


2.1 


1.5 


3.8 


1.3 


2.2 


1.1 


1.9 


1.0 


1.4 


0.8 


0.9 


0.8 


qgsc 


















1.6 


1.0 


1.2 


0.7 


0.8 


0.7 


qgsp 


















1.9 


1.0 


1.2 


0.6 


0.9 


0.8 


ftfp 






















0.8 


0.7 


1.1 


1.0 


mars 


1.3 


1.2 


1.2 


0.9 


0.9 


1.0 


0.7 


0.6 


1.0 


1.1 


0.8 


0.6 


0.8 


1.0 


rpg 


1.4 


0.6 


2.7 


0.8 


1.2 


0.7 


1.8 


0.9 


1.1 


0.6 


0.7 


0.4 


0.9 


0.7 



Generator 

bert 
Ihep 
qgsb 
qgsc 

qgsp 

ftfp 
mars 

rpg 



Be 3 GeV 

1.0 1.1 
0.5 1.6 
1.7 4.8 



1.1 1.2 
0.6 1.7 



Ta 3 GeV 

0.5 0.6 
0.7 2.9 
1.1 1.7 



0.7 1.0 
0.7 2.8 



Be 5 GeV 

1.2 1.4 
0.8 1.9 
1.2 4.0 



Normalization factors, incident tt + 

Ta 5 GeV Be 8 GeV Ta 8 GeV 

+ 



0.9 0.9 
0.7 1.7 



7T+ 

0.9 



7T 

1.2 



1.0 2.6 
0.9 2.3 



0.6 0.6 
0.8 1.9 



7T ' 7T 

1.6 2.1 

0.9 1.7 

1.2 2.2 

1.1 1.8 

1.2 2.2 

1.2 1.2 

0.8 1.5 

1 



7T + 7T 

1.4 1.9 



0.6 
0.8 



1.2 
1.7 



0.7 1.3 

0.7 1.4 

0.7 1.1 

0.7 0.9 

0.5 0.9 



Be 1 2 GeV 

+ 



TT 

1.2 
0.6 



7T 

1.6 
1.0 



0.6 0.8 

0.6 0.7 

0.6 0.8 

0.7 1.0 

0.8 0.7 

0.5 0.9 



